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Abstract 

We compute the total top-quark pair production cross section at the Tevatron and 
LHC based on approximate NNLO results, and on the summation of threshold 
logarithms and Coulomb enhancements to all orders with next-to-next-to-leading 
logarithmic (NNLL) accuracy, including bound-state effects. We find 

att-(Tevatron) = {7 .22tU\-oil) 

aaiLEC, = 7TeV) = (162.6t^:^t}i^) pb 



for mt = 173.3 GeV. The implementation of joint soft and Coulomb resummation, 
its ambiguities, and the present theoretical uncertainty are discussed in detail. We 
further obtain new approximate results at N^LO. 



1 Introduction 



Fifteen years after the discovery of the top quark, the experimental uncertainty of the 
measurement of the top-pair production cross section at the Tevatron has dropped below 
ten percent P-S]. At the LHC, the precision of the measurements quickly approaches 
that reached at the Tevatron OE]. Therefore precision studies of the heaviest quark 
known today can be expected in the near future. An interesting use of the total cross 
section measurement is the extraction of the top-quark mass [7-9j in a cleaner (albeit less 
precise) way than the direct measurements. Furthermore, the total cross section provides 
constraints on new-physics models that try to explain the forward-backward asymmetry 
anomaly observed at CDF [lOj. More generally, its experimental reconstruction involves 
several aspects, such as missing energy and 6-quark tagging, relevant to the search for 
heavy particles in extensions of the Standard Model. 

The precision of experiments motivates improving the accuracy of the theoretical cal- 
culation of the total top-pair production cross section beyond that of the next-to-leading 
order (NLO) result [11], and of the summation of leading-logarithmic (LL) or next-to- 
leading logarithmic (NLL) soft-gluon corrections [T214T7] . Work on the complete next-to- 
next-to-leading order (NNLO) corrections is in progress and several ingredients are already 
known [IBIET] . Meanwhile, progress in the understanding of the infrared structure of mas- 
sive two-loop scattering amplitudes in QCD [2SH33] has provided the prerequisites to com- 
pute soft-gluon effects at next-to-next-to-leading logarithmic (NNLL) accuracy |3l] - H2] . 

In this paper we are concerned with resumming higher-order corrections to the partonic 
cross sections, (jpp/, for the processes 



where the top quarks have a small relative velocity /3 = y 1 In this limit corrections 

due to the radiation of soft gluons and the exchange of virtual Coulomb gluons are both 
enhanced, resulting in terms of the form as log^'^ /3 and as/ (3 at every order of the coupling 
expansion. Near the partonic threshold, /3 ~ 0, these corrections become large and should 
be resummed to all orders. Instead of organizing radiative corrections by the power of 
the strong coupling constant it is then appropriate to count both as\nf3 and as/ (3 as 
quantities of order one. Exploiting the exponentiation of double logarithmic terms leads 
to a parametric representation of the expansion of the cross section in the form 



(1.1) 



that become large in the partonic threshold limit. 



s = {ki + fca)^ 4m^ 



(1.2) 




JJZl^j lnf3go{as In /3) + ^i(a^ In (3) + asg2{as In (3) + . . . 




(LL) (NLL) (NNLL) 

X {1 (LL,NLL); a„ f3 (NNLL); a^, as/3, {3^ (N^LL); . . .} . 



(1.3) 
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We remark that the structure of the expansion is significantly more comphcated than for 
Drell-Yan or Higgs production. The presence of Coulomb enhancements, (cts//?)'^, not 
only calls for a joint summation of soft-gluon and Coulomb effects, but also implies that 
one must address corrections from soft emission at NNLL [36l|39], which count as 

a power-suppressed effect in Drell-Yan and Higgs production. In general, increasing the 
order in (11. 3p requires higher-loop corrections to anomalous dimensions to sum logarithms, 
and the inclusion of power corrections of increasingly higher order. 

In order to obtain the total hadronic cross section, the partonic cross section is convo- 
luted with the parton distribution functions. Both at Tevatron and LHC, the top anti-top 
invariant-mass distribution peaks at about 400 GeV, which corresponds (in the absence of 
radiation) to /3 ~ 0.5. The convolution of the partonic cross section with the parton lumi- 
nosity is therefore dominated by the region /3 > 0.3, where the threshold approximation 
is not valid (see e.g Nevertheless, one often finds that the threshold expansion pro- 

vides a reasonable approximation even outside its domain of validity. For instance, in the 
absence of the exact NLO calculation one would do much better including the threshold- 
enhanced NLO terms rather than staying with the LO approximation. In the same spirit, 
we consider NNLL threshold resummation, which includes the singular terms at threshold 
at NNLO as a subset, a useful step towards the full NNLO result. Even when the latter 
is finally known, the resummed result contains a set of higher-order terms that may shed 
light on the fixed-order expansion of the total cross section. 

The resummation of threshold logarithms in top-pair production has been convention- 
ally performed in Mellin- moment space [ISlIll], and numerical results for resummed correc- 
tions [inillSlllG] at NLL accuracy, or partial higher-order results obtained from expansions 
of the resummed cross section [HIIIT], have been known for some time. In these works the 
Coulomb corrections were included in fixed order, which corresponds to setting k = or 
1 in (11.31) . An approximate NNLO correction based on the expansion of the NNLL result 
up to order has been obtained in [36], while the effect of Coulomb corrections in the 
threshold region has been considered in ^SlBQ] . A combined resummation of soft and 
Coulomb corrections has not yet been performed, and will be done in this paper at NNLL 
accuracy. We employ the momentum-space approach to threshold resummation [5UH5^ . 
and the next-to-leading order Coulomb Green function obtained in the context of top-quark 
production in electron-positron collisions [53l[5l]- We include all corrections at NNLL ac- 
curacy defined by (II. 3p . apart from a class of higher-order non-relativistic (as opposed to 
soft) logarithms, whose summation would require an extension of results of [551156] . The 
fixed-order logarithmic correction at O(a^) was computed in [26] and is included in our 
results. The numerical impact of the resummation of these logarithms is expected to be in 
the per-mille range and their calculation is left for future work. 

It should be mentioned that beyond the higher-order QCD effects discussed in our work, 
electroweak corrections and finite- width effects can become relevant at the level of precision 
aimed at by our NNLL resummation. The former have been computed in [53,158] and are 
of the order of —2% for the total cross section at the LHC, and smaller at Tevatron, while 
the effect on distributions can be larger. Finite-width effects on the total cross section are 
expected to be of the order Tt/rrit < 1%, as has been confirmed by recent calculations of 
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the NLO QCD corrections to the full bbW+W- final state [MIEQ]. 

Recently higher-order expansions and NNLL resummations in different kinematical 
variables, such as the invariant mass distribution and so called one-particle inclusive kine- 
matics [371[3Hlllllll2] , became available. In these calculations different logarithms than the 
ones in the threshold expansion i \1.3\\ are resummed. With respect to the total cross section, 
this amounts to including a particular set of higher-order terms in /3, while excluding some 
of the more singular terms. Since the total cross section receives important contributions 
from outside the threshold region, the two resummation procedures are complementary. 
The numerical results of these studies appear to differ beyond the quoted uncertainties, 
which calls for particular attention to a realistic estimate of the remaining uncertainty 
inherent in the threshold expansion. 

The paper is organized as follows. In Section [2] we review the combined soft and 
Coulomb resummation developed in [3^139] and provide the ingredients required to achieve 
NNLL resummation for top-quark pair production. We expand the NNLL resummed cross 
section to order and present a new result of exact terms beyond NNLO in Section [31 
In Section H] we discuss and define our implementation of resummation and the method to 
estimate the uncertainties resulting from the threshold approximation. Section [5] contains 
our final results for the top-quark pair production cross section at Tevatron and LHC for 
a range of top quark masses. Details of the evaluation of the NLO Coulomb correction, 
including the bound-state contribution, the complete expressions for the expansion up to 
order a^, and a discussion of our error estimate are contained in the appendices. 

2 Basics of resummation 

2.1 Summary of the resummation formula 

Resummation of threshold logarithms is based on the factorization of the partonic cross 
section in the threshold region into several functions [l3l[61] receiving contributions from 
different kinematical regions (soft, hard, coUinear...). The independence of the cross sec- 
tion on artificial scales introduced in this factorization implies evolution equations, whose 
solutions resum threshold logarithms to all orders. For the case of top pairs (and other 
heavy coloured particles) the situation is more complicated, since in the threshold region 
soft gluons and the non-relativistic top quarks have comparable energies of order m^/?^, 
so that energy is modified by a relevant amount due to soft radiation. This is in contrast 
to the case usually considered in deriving factorization formulas at partonic thresholds, 
where all particles are taken to be highly energetic compared to the scale of the soft ra- 
diation. As a result, the factorization formula for the hard-scattering total cross section 
for the partonic subprocesses (11.11) near the partonic threshold (II. 2p factorizes into three 
contributions [39], with a potential function J in addition to the hard and soft functions 
H and W: 
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Here E = \/§ — 2mt is the energy relative to the production threshold. The sum is over the 
colour representations of the final state top-pair system, i.e. the colour-singlet and octet 
states. The hard functions H can be obtained from the colour-separated hard-scattering 
amplitudes for the process pp' — )■ tt in the threshold limit, as described in |39j. The 
soft function is defined as a vacuum expectation value of soft Wilson lines, and can be 
reduced to one for a fictitious 2-^1 scattering processes with a single final-state particle in 
the colour representation R, provided an appropriate s-channel colour basis is chosen [31]. 
Such a basis has been adopted in (12. ip . The potential function sums Coulomb gluon 
exchange related to the attractive or repulsive Coulomb force in the colour singlet and 
octet channels, respectively. 

In the momentum-space approach [501 - I52] resummation is performed by evaluating the 
hard function H{fi) at a hard matching scale fih ~ 2mf and the soft function W{uj,fi) at 
a soft matching scale of the order of /ig ~ mtP'^. Both are then evolved to an intermediate 
factorization scale using evolution equations. For the case of heavy pair production at 
threshold the relevant evolution equations and the solutions have been given in [MlEn] , so 
we will only summarize the solutions below. 

The convolution of the soft and potential functions in (12. ip accounts for the energy 
loss of the top pair into soft radiation. In our study of NLL resummation for squark anti- 
squark production [39] we observed a non-negligible numerical effect from the combined 
soft and Coulomb resummation, and also found a reduced scale dependence. For top anti- 
top production the Coulomb corrections are smaller, so we expect less sizeable effects for 
the case of top quarks. 

The formula (12. ip has been derived in an effective field theory for the production of 
heavy particles in an S'-wave stat^ and is valid for the leading effective Lagrangian. Correc- 
tions from subleading terms in the effective Lagrangian, or higher-dimensional production 
operators, potentially lead to the corrections 0{(3) (NNLL) and 0{as(3, (N^LL) indi- 
cated in (II. 3p . and would have to be accounted for by including additional hard, potential 
and soft functions in (12. ip . As shown in [36l|39], these subleading effects vanish at NNLL 
accuracy due to rotational invariance, so that the 0{P) corrections are absent and factor- 
ization holds in the simpler form (12.11) . Starting from N^LL order, such power-suppressed 
corrections in /3 should be expected, requiring the introduction of new soft and poten- 
tial functions, as is also indicated by the massive two-loop soft anomalous dimension for 
non-threshold kinematics [3T1[32], which exhibits a more complicated colour structure. 

2.2 Inputs to the resummation formula 

From (II. 3p and the factorization formula (12. ip . it can be seen that the ingredients required 
for NNLL resummation are given by the resummation functions in the exponent [39] [52]. 

^Therefore the formula applies to the qq and gg initiated partonic channels but not to contributions 
from qg, qg initial states, that are suppressed by a factor and therefore formally of N'^LL order. Results 
for the leading threshold-enhanced NNLO contributions for those channels have been given in j7j. 
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the fixed-order one-loop soft [M] and hard [62] function^ and the NLO-resummed poten- 
tial functions [53]. Once the only non-trivial pro cess- dependent quantities, the one-loop 
hard functions, are known, the formalism can be immediately applied to perform an NNLL 
resummation for other processes, such as the production of coloured supersymmetric par- 
ticles. In this subsection we collect the required ingredients. 



2.2.1 Hard functions 



The hard functions are related to the amphtudes of partonic top production processes 
pp' — 7- ti directly at threshold. Their perturbative expansion in the QCD coupling in the 
MS scheme can be written as 



(2.2) 



The leading-order hard functions are related to the partonic Born cross section at threshold 
in the colour channel R, a^J, ^, according to 



.(0) 



27r 



(2.3) 



The NLO coefficients h^^, at the hard matching scale fih can be obtained by comparing 
the analytical result for the threshold expansion of the NLO top-pair cross section [62] to 
the NLO expansion of the resummed cross section [39]. Subtracting the NLO contribution 
from the Coulomb and soft function, we obtain 
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Ami 



+ 2Cf 
57r2 



— - 10 



(2.4) 



12 



+ 2C, 



^-10 



in the MS scheme with five active quark fiavours. Here R = 1 stands for the colour- singlet 
and i? = 8 for the colour-octet channel, respectively. The colour factors are Ca = Nc = 3, 
Cp = {N^; — l)/{2Nc) = 4/3. While strictly speaking the leading term in the threshold 
expansion of the Born cross section appears in the hard function (12. 3p . in our numerical 
results below we use the exact expression for the Born cross section, that can be found e.g. 
in 



^The 0{as) constant contribution appearing at NNLL in (|1.3p has been included in some earlier NLL 
predictions |16| in an approximate way. 
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The choice /i/i ~ 2mt ehminates the logarithmic terms in the hard function at the 
matching scale, but the soft and Coulomb functions naturally live at parametrically smaller 
scales, of the order of and mj/3, respectively. The solution to the renormalization 

group equation that sums logarithms of ln(/i/i//i) is given by 

= exp[4^(/i,,/i) - 2ar(/x,,/x)] , (2.5) 

V P/i / 

with the functions S, af and ar defined as [52] 

"^^'^^ , r^usp(«s) + rlp(a.) da'. 



S{jj,h,iJ,) = - I da 



ar(/ih,/i) = - da 

J a.,(uh) 



2/3(«.) 

ar(/i.,/i) = - / ^^s^j^. (2.6) 

The labels r, r' denote the colour representation of the initial-state partons p, p', while the 
label % refers to the colour of the initial-state partons and the representation R of the tt pair. 
Explicit results for the MS /3-function, ap and the Sudakov exponent S up to N^LL order 
can be found in [52]. The relevant anomalous dimension coefficients as needed for NNLL 
accuracy are collected in [Mj. Starting from the two- loop order the evolution equation of 
the hard function is modified because of additional IR divergences, that are related to UV 
divergences of the potential function due to insertions of non-Coulomb potentials. In this 
work we will not consider these contributions (see end of section [2.2.31 below). 

2.2.2 Soft functions 

The resummed soft function in the momentum-space formalism [501452] can be written as 

1 f p-27-E'? 

Wf'^\u,y) = exp[-45(/i„/i) + 2a^i/i„/i)] Sf (9,,/i.)- - 6(00)— — . (2.7) 



f^sj ' 'r(277) ■ 



The auxiliary variable 77 is set to rj = 2ar(ps)Ai) after performing the derivatives. The 
function a^,^ is defined, analogously to in (12. 6p . in terms of the anomalous dimension 
■j^r^i, first obtained up to the two-loop level in [31]. The function sf (p) is the Laplace 
transform of the soft function. For resummation at NNLL accuracy we require the NLO 
soft function, which reads [311155] 

\o{al), (2.8) 



(a +a.') (p' + y ) -2Ch(p-2) 



where Cr is the quadratic Casimir operator in the colour representation r of the initial 
state parton p. As mentioned above, the soft matching scale p^ should be chosen of the 
order m^/?^ in order to resum large logarithms in the soft function. The precise choice we 
adopt in our numerical analysis will be discussed in Section |H 
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2.2.3 The NLO Coulomb function 



The potential function is related to the imaginary part of the zero-distance Coulomb 
Green function of the Schrodinger operator — V^/mt — {—DR)as/r. The coefficients of 
the Coulomb potential in the colour- singlet and colour-octet configuration are given by 



2N, 



c 



2N, 



(2.9) 



c 



For NNLL resummation we require the Coulomb Green function up to next-to-leading 
order, which includes one insertion of the NLO Coulomb potential 



5V{p,q) 



ai — /3o In ■ 



(2.10) 



Y^A — ^^iTf is the one-loop beta-function coefficient, and ai = yC^ — "^riiTf. 



where /3o 

The insertion of 5V(j),q) yields a factor x ag/P, which according to (11. 3p produces a 
NNLL correction. The potential function up to next-to-leading order can be written as 



Jr{E) = 2 Im Gg)^(0, 0; E) Kc{E) + G^^^^(0, 0; E) + . . . 



(1) 



(2.11) 



where G^^^ is the solution to the Schrodinger equation with the leading Coulomb potential, 

resumming all (a,,//?)" corrections, while G^^'^j^ sums x (a,,//?)" corrections by solving 

perturbatively the Schrodinger equation with one insertion of 5V{p^q) [53j. The explicit 
expressions are given in Appendix lAl 

The NLO potential function is independent of the scale /i, though a residual scale 
dependence, formally of higher order, remains. We use a /3-dependent scale [39] . 



Ixc = max{2mi/3, Cprntasific)} 



(2.12) 



in the numerical evaluation. With this choice logarithms of (3 in the potential function 
arise only from the scale of the strong coupling. The lower limit in f l2.12p is chosen to 
coincide with the Bohr scale of the lowest-lying bound state in the colour singlet channel, 
which is indeed the relevant scale for /3 — )■ in this case, as can be seen from the argument 
of the logarithm in (lA.SP below. There are no bound states in the repulsive colour-octet 
channel, but since the Coulomb effects are small in this case, the precise choice of scale is 
not important. 

Further NNLL efi'ects from the non-relativistic dynamics arise from NNLO terms in the 
heavy-quark potential not related to the Coulomb potential. 



NNLO 



[p, q) 



7ias{fi'^)\q\ 
4m 



Dr 

2 



Ca 



+ 



m 



2 ^spin 



(2.13) 



where Vspm = and —2/3 for a tt pair in a spin-singlet and spin-triplet state, respectively. 
The evaluation of an insertion of the non-Coulomb potential results in UV divergences that 
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cancel against IR divergences in the two-loop hard function, as mentioned in Section 12.2. 1[ 
After this cancellation a logarithm of f3 remains at NNLO, which counts as NNLL. The 
divergence also produces an additional contribution to the anomalous dimension •yY of 
the hard function starting from the two-loop order. Since as mentioned above, we do 
not consider this additional contribution in the NNLL resummed hard function, the non- 
Coulomb corrections will be included only at fixed order in the number of non-relativistic 
logarithms. The corresponding NNLO correction for a final state in an arbitrary colour 
representation is given by [36] 

ax\nc = aP alifx) In (3 [-2Dl (1 + t;spin) + DrCa] • (2.14) 

Here a^'' is the Born cross section in the spin and colour channel X. For top quarks the 
Born cross section in the qq initiated channel is a pure colour-octet spin-triplet, whereas in 
gluon-gluon fusion the ti state is spin-singlet but colour-octet or singlet. In our numerical 
results, the corrections ( 12.14P are added to the potential function Jr{E) in (12. lip through 
the factoiil 

A,c(^) = 1 + «'(/ic) ln/3 [-2DI (1 + v,^,^) + DrCa] e{E) (2.15) 

multiplying the leading Coulomb function G^^j:j(0, 0; E)^ with (3 = {E/mtY^'^. After convo- 
lution with the soft function this procedure yields the correct In (3 and In^ /3 terms in 
the threshold expansion. The factorized form of the logarithmic non-Coulomb correction 
can be deduced from the results of [53] (given explicitly in the appendix of [SB]) and [M] . 
and produces correctly the series a^log/3 x {as/PY of terms associated with the first 
non-Coulomb logarithm and any number k of Coulomb exchanges. 

The non-relativistic anomalous dimension produces a series x (a^log/?)"" x (as/f3)^, 
of which (I2.15P is only the first term (n = 1, any k), that should also be summed for 
complete NNLL accuracy. This requires the generalization of results for top-quark pair 
production in e~^e~ collisions (see, e.g., [55l[56]) to the production of a colour-octet final 
state. The numerical effect of the leading non-Coulomb terms is an enhancement of the 
top production cross section of about 0.5% at Tevatron and LHC. The higher-order terms 
should therefore be very small, and we leave this resummation of non-relativistic logarithms 
to future work. 

2.3 Calculation of the convolutions 

The total hadronic cross section for the production of a tt + X final state in collisions of 
hadrons A^^i 2 with centre-of-mass energy s is obtained from 

(TNiN2-^ttx{s) = ^ / dTLpp>{T,fif)app'{sT,fIf) , (2.16) 
p,p'=q,q,9 -^4"^?/^ 

•^For simplicity we implement the non-Coulomb correction factor only in the continuum, and not for 
the bound states, since the uncorrected bound-state effect is already very small; hence the theta-function 
in 1135]). 
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where the parton luminosity is defined in terms of the parton distributions functions (PDFs) 



Lpp'{T,n)= / dXidx2 S{XiX2 - t) fp/N^{xi, H)fp> /N2{X2, IJ') 



(2.17) 



We briefiy discuss some technical issues encountered when performing the convolution of 
the resummed soft and potential functions in the partonic cross section (12. ip . and the 
subsequent convolution with the parton luminosity. 

We convolute the NLO Coulomb Green function (12. lip with the resummed soft function, 
including the (pseudo) bound-state contributions from energies below the nominal top anti- 
top threshold. We note that for integrated quantities, such as the total cross section, the top 
width can be set to zero in the Coulomb function. Thus, for the continuum contributions, 
the integral in (12. ip has to be evaluated for < u < 2E. The bound-state contributions 
are located at values Un > ^E, and we discuss them in Appendix [B] Here we focus on the 
continuum. For the relevant hierarchy of scales, fig < f^f, the anomalous dimension r] is 
negative and the expression ( 12. 7p should be understood in the distributional sense [H]. 

Let us first discuss the simpler case where the non-Coulomb corrections to the potential 
function are neglected, i.e. Anc = 1 in (12. lip . The convolution of the soft and potential 
functions in the factorization formula (12. ip is analytically continued to negative values of 
T] according to 



2E 



dujjR^E - — 

















Jo w - 



+ 



ME - |) 



ME) ME)E 



27] 



2r] + l 



ME) + ^ME) 



2-n 



(2.18) 



The derivatives of the potential function appearing in (I2.18P have been evaluated by ex- 
pressing the hypergeometric function in terms of harmonic sums using methods from [55] , 
as described in Appendix [XI The double subtraction in (I2.18P is sufficient to render the u 
integral finite for rj > —1. However, for rj < —1/2 the convolution of the first term in the 
second-line square brackets with the parton luminosity in (I2.16P requires regularization at 
the partonic threshold r — )■ Amf/s = tq, if the potential function is non- vanishing at E = 0, 
as is the case for the attractive colour- singlet Coulomb potential. To see this, note that in 
the region r ~ tq we can approximate E = y/sr — 2mt ~ ^y~{t — tq). With L 



pp. yT, fif) 



and J{E) approaching constants as r — )• tq, the r-integrand behaves as (r — tq)^^. 



Thus 



the r-integral for the colour- singlet cross section is analytically continued to rj < —1/2 
according to 



^ dTLpp,{T)J,iE) (^^^ ' = —L-Lpp,{To)UO){l 



s — Amf 
2mtfis 



2\ 2»7 



1 

+ / dr 

'TO 



Lpp,{T)Ji{E) 



2E 

/is 



2ri 



- Lpp,{To)Ji{0) 



ST — Arrii 
2mtfis 



2\ '^v' 



. (2.19) 



9 



The pole at = — ^ in the first hne cancels with the overall gamma function in the 
resummed soft function (12. 7p . so that the cross section (12.161) is now well-defined for > — 1. 
In order to avoid large numerical cancellations, it is practical to split the integration interval 
into two regions and perform the analytical continuation (I2.19P only in a small region 
r G [tq, To + At]. The potential function in the singlet channel a.t E = 0, Ji(0), is given 
in flA.5|) . The derivative of the singlet potential function at E = vanishes only at leading 
order so that further subtractions would be required when fig and /i/ are such that t] < —1. 

Alternatively, the left-hand side of (I2.18P can be analytically continued to negative 1] by 
performing a partial integration with respect to u. After one partial integration we obtain 



/ dooME--) 
Jo ^ 



1 (u 



2r, 



1 

'2r] 



2E 



duj 



dJR{E - 
du 



2t] 



(2.20) 



The integral in the second line is well defined for > —1/2. By applying a second partial 
integration, it can be continued to values > — 1, as in fl2.18p . generating boundary terms 
with derivatives of the Coulomb function. The boundary terms at a; = are zero in 
the analytic continuation from positive t]. For those at a; = 2E only the singlet channel 
contributes, since the octet Coulomb function Js{0) and its derivatives at -E = vanish. 
The convolution of the boundary terms with the parton luminosity can be treated in 
the same way as the terms in the second line of (I2.18p . leading to expressions similar 
to (I2.19p . Note that for the LO Coulomb function only Ji(0) is non-zero, whereas all 
higher derivatives vanish. 

Including the non-Coulomb correction (12.150 in the potential function (12.110 makes 
the analytic continuation more complicated due to the additional term J^\e) \og{E/mt), 
where J^\e) = Im [G^^^(0, 0; i?)] . The octet case does not pose a problem, since the 
octet Coulomb function and its derivatives vanish faster than any power as — )■ 0. But for 
the sing let ji°\E)\n{E/mt) diver ges in this limit, and we have to introduce an additional 
subtraction, 



2E 



nit 



duJi'"{E--)\n 



2E 



+ I doj[jf\E- 



2r? 



Jf(0) 



2E 



du In 



E-J 
mt 



4»)(0)),n(^ 

/ \ rrit 



2n 



2rj 



(2.21) 



The first integral on the right-hand side of (12.210 can be calculated analytically, whereas 
the second integral can be continued to lower values of t] either using partial integration, 
as in or as in fITTH]) . with the replacement Ji(E) ^ (jf^(^) - jf ^(0)) ln(E/m). 

The convolution of the luminosity function Lpp/(r) with the second integral on the 
right-hand side of (12.210 is finite, once the u integration has been continued to negative 
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values of rj. On the contrary, the first integral, 

,<«.,o) r\. (^f . L (A) _ _ , , 

(2.22) 

contains a non-integrable singularity at = for rj < —1/2. This can be subtracted and 
analytically continued to rj > —1, similarly to (I2.19p . 

3 Expansion to 

The expression for the resummed cross section (12.11) can be expanded to order a" in 
the strong coupling, providing an approximation to the full C(a^) QCD calculation. The 
expansion of the NLL and NNLL resummed corrections up to order a J and a^, respectively, 
has been given in [36|[39]. Here we provide the corresponding expansion of the NNLL result 
up to order a^. In this way, we generate the first terms in the threshold expansion of the 
third-order fixed-order cross section, whose size should be indicative of the quality of the 
fixed-order expansion. If there are no cancellations, integrating the threshold expansion 
gives a rough estimate of the correction that may be expected on top of the exact NNLO 
result, once it is known. 

To generate the threshold expansion from the resummation formula, the scales have to 
be chosen as fig = ksmt/3'^, fic = kcmt/3, fih = khrut, where the kx are numbers of order 
one. Those terms in the threshold expansion that are completely included at NNLL are 
independent of the kx- The expansion of the resummation formula also includes terms 
that are beyond the NNLL accuracy, which may depend on the arbitrary constants kx, 
and are presently not completely known. Adopting the notation 

V,H(A./) = *:!,«{l + EE(^)"4"5,ln'"(^) } (3^1) 

for the fixed-order expansion in the strong coupling of the partonic cross section in colour 
channel R, the result for the different production channels can be written as 

/^^Jgj =12945.4 In^ /3 - 37369.1 In^ /3 + 27721.4 In^ /3 + 41839.4 In^ /3 

+ ^ [-2994.51 In^ /3 + 2804.73 In^ /3 + 3862.46 In^ (3 - 6528.61 ln/3] 

+ [153.93 In^ /3 + 122.866 In /3 - 144.996] + f^^-'^^ , (3.2a) 
=147456. In^ (3 - 59065.6 In^ /3 - 286099. In^ (3 + 349463. In^ (3 
+ i [121278. In^ (3 + 103557. In^ f3 - 164944. In^ (3 + 56418.5 In (3] 

+ [22166. ln2/3 + 39012.1 In /3- 2876.61] + /^gfij , (3.2b) 
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=147456. In^ /? - 169658. In^ /3 - 140834. In^ /3 + 524210. In^ /3 

+ ^ [-15159.7 In^ /3 - 5364.82 In^ (3 + 19598.9 In^ (3 - 17054.7 In f3] 

+ [346.343 In^ /3 + 522.978 In /3 - 71.7884] + /^J^gJ . (3.2c) 

The terms given explicitly in (13. 2p are those N^LO terms, which are predicted correctly by 
the NNLL approximation, and which belong to the NNLL terms according to (11.31) . The 
remainder functions fi^'^\ containing N^LL and higher-order terms, are collected in Ap- 
pendix [Cl together with the coefficients of the \n{fif/mt) terms. There we also summarize 
the O(a^) terms from |36J and the /cx-dependent constants not given in [36], which are 
generated from the expansion of the NNLL resummation formula, but are formally beyond 
NNLL. 

We note that while the Coulomb terms are generically 0{{as/ P)"'), there is no l//3^ 
term at third order. Instead a bound-state contribution is present, that in the fixed-order 
expansion would appear as a S{/3) term, which is not included in (13. 2 p but can be deduced 
from the expressions given in Appendix [Bl The absence of a 1/(3^ term implies that the 
convolution of the partonic cross section with the parton densities is still well-defined, since 
dcTpp'/d/S (X In^ l3 for small f3. This ceases to be true at N'^LO, where a 1//?"^ term is present 
in fi^'^\ At this order the fixed-order expansion of heavy particle pair production breaks 
down, and the summation of Coulomb terms is mandatory even for the total cross section, 
which is not necessarily threshold-dominated after the resummation has been implemented. 

In order to judge the magnitude of the (9(ag)-terms, we plot the corrections to the 
partonic cross section in Figures [1] and O We consider three approximations: 

• ANNLOapp, which consists of the singular terms of the C(a^) correction as /3 — )■ 
[36]. 

• N^LOa, which adds to the above the 0{al) terms from (13. 2p . including the remainder 
functions. 

• N^LOb, which adds to ANNLOapp the 0{al) terms from (13. 2 p with the remainder 
functions set to zero, i.e. the NNLL terms only. 

The difference between the first and the other two approximations provides an estimate 
of the importance of the third-order terms, and hence the convergence of the fixed-order 
expansion. The difference between the two N^LO approximations represents an ambiguity 
in our estimate of the third-order termsEI 

In the figures we plot the integrand of the convolution in the formula for the hadronic 
cross section (I2.16p . i.e. the product of the partonic cross section and the parton luminosity, 

"^For this estimate we choose kg = kc = 1 and kh — 2, which differs from our canonical choice 
ks = kc = kh = 2 discussed below. The reason for this is that the remainder functions turn out to be 
small for the canonical choice, so the difference between N'^LOa and N'^LOb would not provide a measure 
of the "natural" size of the sub-leading singular terms at N'^LO. 
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Figure 1: Partonic higher-order cross sections for qq — t- tt (top) and gg — t- tt (below) at the 
Tevatron, multiplied with the parton luminosities. Black (solid): ANNLOapp, blue (dashed): 
N^LOa (all terms in for ks = kc = 1, kh = 2), red (dot-dashed): N^LOb- 



function of /3, including the Jacobian dr/dfS: 

;Lpp'{l3,iJ,f)Aapp>^tt{l3,iJ,f) . (3.3) 



d/3 ~s(l-/32)2- 

We show the partonic cross sections at the Tevatron (Figure [1]) and LHC {^/s = 7TeV, 
Figure |2]), using the MSTW2008NNLO PDF set [64J and = rut. The third-order cor- 
rections are not neghgible. The strong increase at small /3 should be expected, since for 
such small values of /3 the perturbative expansion breaks down and resummation should be 
performed. However, we note that for the gluon channel the magnitude of the third-order 
terms can be comparable to the second-order ones up to /3 ~ 0.6 (N^LOa). Unless there 
are cancellations with higher-order terms in the (3 expansion, which need not be small at 
P ~ 0.5, this may indicate a poor convergence of the fixed-order perturbative expansion. 
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Figure 2: Partonic higher-order cross sections for qq — )■ tt (top) and gg — )• tt (below) at the 
LHC (a/s = 7TeV), multiplied with the parton luminosities. Black (solid): ANNLOapp, blue 
(dashed): N^LOa (all terms in for ks = kc = 1, kh = 2), red (dot-dashed): N^LOb- 

The area under the curves in the figures gives directly the correction to the total 
hadronic cross section that needs to be added to the NLO prediction. Adopting the 
top mass value mj = 173.3 GeV, we obtain an additional 0.60 pb for ANNLOapp, 0.93 pb 
(0.71 pb) for N^LOb (N^LOa) at the Tevatron and 12.1 pb for ANNLOapp, 17.1 pb (7.8 pb) 
for N^LOb (N^LOa) at the LHC with ^/s = 7TeV. The third-order term alone may 
therefore amount to up to 4% of the cross section at Tevatron and LHC. This might be 
considered as an estimate of the intrinsic uncertainty of the fixed-order prediction after the 
exact NNLO result is known. Below we shall see that resumming the threshold expansion 
to all orders in the NNLL approximation leads to a smaller effect from the sum of all 
higher-order terms. 



14 



4 Implementation of resummation and uncertainties 



In the implementation of the resummation formahsm one encounters several ambiguities: 
kinematical ambiguities, the precise matching to the fixed-order calculation, the choice of 
the scales fih and fis appearing in the momentum-space approach to threshold resummation, 
and the scale /ic in the Coulomb function. We here discuss the choices made and the 
procedures to estimate the resulting theoretical uncertainties. 



Kinematic ambiguities. If the total hadronic cross section fl2.16p is computed by con- 
voluting the factorized partonic cross section (12. ip with the parton luminosity for arbitrary 
values of r, the application of the factorization formula is extended outside its domain of 
validity s ~ 4mf. As a result, there are kinematic ambiguities, since expressions that agree 
at partonic threshold can show numerical differences for /3 — )■ 1. One such ambiguity is 
related to the expression used for the hard function. As already mentioned in Section !^. 2. 
we use the exact Born cross section, instead of the threshold limit, in the determination of 
the leading hard function (12. 3p . since this choice is observed to lead to a better agreement 
of the threshold approximation with the exact NLO result. Another ambiguity arises from 
the argument of the potential function J^j in (12. ip . The derivation of the factorization for- 
mula [39] results in the dependence on the exact energy relative to the nominal production 
threshold, E = y/S — 2mt = 2mj x (1/ ^yl — /S"^ — 1). However, the threshold expansion 
is customarily formulated as an expansion in /3. In the threshold limit the two choices 
are equivalent, since E ^ Tnt/S"^, and hence one may replace E by m^/?^ in the potential 
function. We have implemented both choices but use the second one as our default, since 
in this case the expansion parameter E/rrit ^ (3"^ is always smaller than 1, whereas it can 
grow to large values when the exact expression for E is used. The difference between the 
two implementations is used as a numerical estimate of the ambiguities in the resummation 
procedure. 



Matching to the exact NLO result. Since the total production cross section is not 
numerically dominated by the threshold region, where the factorization formula (12. ip pro- 
vides an accurate description of the cross section, we match the NNLL calculation to the 
exactly known NLO result. We consider two options. In the first, the NNLL cross section 
is expanded to NLO and subtracted from the resummed result. The resulting higher-order 
corrections are added to the full NLO result: 



"-^pp', matched, 1 V^J 



a: 



NNLL/ 



PP 



[s) - a,^ 



NNLL(l) , 



PP 



(4.1) 



This matching prescription will be denoted by NNLLi in the following. In this way, there 
is no double counting of 0(as) corrections. The expansion of the NNLL-resummed cross 
section to 0(as), ci^^^^^^\ coincides with the approximate NLO cross section given e.g. in 



(B.6) of For the NLO cross section we have implemented both the analytical result 
and the parameterization given in [OS]- For the former, we used the program provided in 
for numeric evaluation. 
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The matched cross section (14.11) contains a constant term at NNLO, i.e. a /3-independent 
correction multiplying the Born cross section. This includes the product of one-loop hard 
functions with the constant terms in the one- loop soft function (12. Sp . as well as terms 
related to ambiguities in the choice of the various scales. Since the two-loop hard and soft 
functions have not been computed yet, the constant NNLO term is not known completely, 
so the inclusion of a partial result is a matter of choice. We therefore consider a second 
matching option, denoted by NNLL2, that matches to the NLO result by subtracting the 
expansion of the NNLL corrections to Oia^) and adding back the NNLOapp corrections: |36] 
with the unknown constant set to zercu 



-NNLL (Q\ 
'^pp'matched,2l''J 



^NNLL(g)_^NNLL(2)(^^ 



PP 



PP 



, -NLO/s\ , -NNLO 



(4.2) 



The expansion of the NNLL correction to (9(a^), a, 



NNLL(2) 
pp' 



is given in Appendix O The 



numerical difference between the two implementations (14. ip and (14. 2 p can be considered as 
an estimate of the unknown constant term at NNLO. We discuss the issue of estimating 
the unknown constant further below. 

For comparison, in the following we also present results for resummed cross sections 
with NLL accuracy. In addition to making the appropriate truncations for the resummation 
functions (12. 6p . in this approximation only the leading order terms in the soft (12. 8p and hard 
function (12. 2 p are used, and the higher-order Green function G^^^j^ and the non-Coulomb 
corrections Anc are set to zero in the potential function (12. lip . In contrast to the NNLL 
predictions, bound-state corrections will not be included in our NLL results presented in 
this paper. Here the two options to implement the matching are given by NLLi, defined 
by 

(4.3) 



cr, 



^NLL I 

pp', matched, 1 > 



-lNLLz-n 
^pp' (SJ 



-NLL(0)/-x 



and NLL2, defined by 



cr, 



NLL 



PP 



.matched, 2 \ 



^NLL(-)_ -NLL(1)(^I 



PP 



PP 



I -NLO /AN 
+ (Tpp, (S) 



(4.4) 



Because of our choice of the hard function, the expansion of the NLL correction to order 
coincides with the exact LO cross section 6"il9(s) 



a^LL(O) 
"s; "pp' ; 



PP 



SO that 0"^'^,matched,l("^) 
NLL(l) 



(T^V^(s). The expansion of the NLL correction to 0{al), denoted by o-^pT"'"', is given 
in [22] ■ The NLL2 matching is identical to the prescription used for squark-antisquark 
production in [39]. Since the full NLO corrections to top-pair production are known, this 
is the preferred implementation for phenomenological results, whereas the NLLi option is 
used to estimate ambiguities in the resummation procedure at NLL. 



Unknown constant at (9(a;^). As mentioned in the previous paragraph, the constant 
term in the threshold approximation is unknown at (9(a^). As a default, we do not include 

^For consistency the constant terms are set to zero also in the coefficients fx''^^ ^nd fx'^^ of ttic 
factorization-scale dependent terms, even though they are known in this case, see Appendix lC.il 
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such a constant in the results, unless it is generated by the resummation, as in the NNLLi- 
option discussed below (14. ip . However, we include an estimate of the order of magnitude 
of the constants in the determination of the uncertainty. Adopting the MSTW08NNLO 
PDFs, fif = rrit and rrit = 173.3 GeV, which is our standard choice, the numerical effect of 
non-vanishing constants for the different partonic channels is given by 



Tevatron : Acx- 
LHC(y^ = 7 TeV) : Ad^ 
LHC(v^ = 14 TeV) : Aa^ 



(2) 



(2) 
tt 



(2) 
tt 



0.345 



1.70 



5.34 



(2) 
QQ 



1000 



a 



(2) 
QQ 



1000 



^(2) 

1000 



^(2) \ 

+ 0.024 +0.008 
' 1000 ' 



'^(2) \ 

+ 4.31 ( + 1.31 

1000 



1000 J 



Pb: 



27.14 



I j 



1000 j 



7.97 



r-(2) 

1000 



^(2) 

1000 



pb, 

pb . 

(4.5) 



Here our conventions are such that the correction to the partonic cross section due to the 
constant reads 



Aa, 



(2) 

pp' ,R 



0", 



(0) 

pp' ,R 



An) 



(2) 

PP',R ' 



(4.6) 



where a^^J, ^ is the (exact) Born cross section for the pp' initial state in the colour channel R. 



pp',R 



(2) 

In order to get a crude estimate of the order of magnitude of the unknown C^p, ^, we 
consider the square of the corresponding one-loop constants. This method is observed to 
give a conservative estimate in cases where the full two-loop corrections are known. For 
top-pair production the one- loop constants have been calculated analytically [62]. In our 
conventions, they are related to the hard functions (12. 4p by 



a 



(1) 

pp',R 



hy{mt)+4{Cr + Cr 



91n^2 - 121n2 + 8 - 



llvr 
"24" 



21 



12Cr [In 2-1] 



Numerically, 



eg) = 14.531 , 



30.586 



*^gg,i 



14.026 , 



hence from \C^]'^'^' 



(Cpp! we obtain the estimates 



I gg 



211.2 



|^(2)est. 

l*-^gg,8 



935.5 



|^(2)cst. 

I^^gg.i 



196.7. 



(4.7) 
(4.8) 

(4.9) 



In our numerical results we estimate the uncertainty due to the unknown constant by 

-^(2)' 



adding and subtracting \C^,'^^'\ in all partonic channels. 



Choice of the soft scale. As discussed in Section |2l the resummation of threshold log- 
arithms in the momentum-space approach amounts to resumming logarithms of the ratios 
and /ifiZ/i/ of the hard, soft and the factorization scale. While an all-order solution 
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to the evolution equations would be independent of these scales, a residual dependence 
remains after truncation of the perturbative series of the resummation functions S and 
a. Furthermore, the potential function depends on the Coulomb scale fic- Our default 
choices for the hard scale, Hh = '2mt, and the Coulomb scale (12.121) have been explained in 
Section [221 

The proper choice of the soft scale requires some discussion. For any given /3, the 
idea of scale separation in the effective theory suggests the /3-dependent scale fig ~ ^tP"^- 
Indeed, this choice has to be made to recover the threshold logarithms in fixed-order in 
as from the expansion of the resummed result. However, the use of a running scale in the 
resummed partonic cross section is problematic, since it leads to an oscillating partonic 
cross section for small values of (3. The problem arises from i] = -^^^^^ log(^) + ... 
becoming increasingly negative, such that the factor l/r(2?7) in (12.71) changes sign whenever 
a pole of the Gamma function is crossed. This effect appears in resummation for all hadron- 
collider processes such as Drell-Yan and Higgs production. In the traditional approach to 
threshold resummation in Mellin space, a similar difficulty appears when performing the 
inverse Mellin transform (see e.g. [13]), since the convolution of the resummed partonic cross 
section with the parton luminosity does not converge. The most widely used prescription 
to perform the convolution was proposed in [13], and amounts to performing the inverse 
Mellin transform at the level of the hadronic cross section. Other prescriptions have been 
suggested to perform the inverse Mellin transform directly for the partonic cross section 
but, in one way or another, employ a cutoff in the convolution with the parton luminosity, 
see e.g. [HI [68]. 

We investigate two methods to determine the soft scale: (1) a fixed soft scale determined 
according to a procedure proposed in [52] and (2) a running scale frozen to a fixed value 
below a certain value of /3. In our final analysis we use the second method by default. 



4.1 Method 1 

In [501452] it has been advocated to choose a fixed soft scale fig- In this approach one does 
not aim at predicting the partonic cross section locally as a function of /3, but argues that 
logarithms in the hadronic cross section for tq = Amf/s — )■ 1 are resummed correctly by 
a fixed soft scale Hs = 2mtks{l — tq), where ks ~ 0.1 due to the steeply falling parton 
luminosity function [52]. The partonic cross section is then predicted only in an average 
sense J§ 

In practice, the fixed soft scale is chosen such that it minimizes the relative fixed-order 
one-loop soft correction to the hadronic cross section. Following this procedure, we vary 
the scale in the PDFs and the soft correction and determine the value ftg that minimizes 
the relative soft corrections: 

= /isT^> / dTLpp,{T,fis) ^^.Q. — . (4.10) 



^See also [42] for a discussion of the differences between this approach and the use of fixed-order 
expansions. 
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Here the fixed-order NLO soft correction, o'^J/gf^f^, is obtained from tlie tliresliold expansion 
of the NLO cross section (see e.g. (D.3) in [39j) by setting the Coulomb correction and the 
hard corrections hi to zero. In the denominator we divide by the leading-order hadronic 
cross section. 

In this approach the value of the soft scale depends on the interplay of coefficients of the 
logarithms and constants in the one-loop soft function fl2.8p . which might be considered 
as going against the philosophy of resumming only large logarithms. Furthermore, the 
reduction of the soft scale by a factor of 10 due to the behaviour of the PDFs compared to 
the naive expectation ~ 2mt(l — tq), has been estimated [52] to be effective for tq > 0.2, 
while for tq — )■ the relevant soft scale is expected to be of the order Hs ^ mt. For top- 
quark production at the Tevatron and the LHC with centre-of-mass energies up to 14 TeV 
the relevant range is tq = 0.03 — 6 x 10~^, so it cannot be expected a priori that the effect 
of the higher-order log /3 terms in the partonic cross section can be properly incorporated 
in an average sense through a fixed /ig. 

As discussed previously, there is a kinematic ambiguity from expressing the one-loop soft 
corrections in terms of \og{E/fis) or of log(mt/3^/ fig). For a process dominated by threshold 
dynamics, the effect of this ambiguity on the soft scale determined through fl4.10p should 
be small. However, for the case of top-quark production, we obtain noticeable differences 
between the soft scales determined using the two parameterizations of the soft function: 



\og{E/fis) : /i« = 52 GeV (Tevatron), 99 GeV (LHC7), 120 GeV (LHC14), 
log(mt/3V/i,) : /i, = 35GeV(Tevatron), 58GeV(LHC7), 65 GeV (LHC14). 



This difference, up to a factor of almost two, is a result of the fact that the total cross 
section is dominated by values (3 > 0.3, and is another hint that the use of a fixed soft 
scale is somewhat problematic for top quark production^ 

We nevertheless present results obtained using a fixed soft scale, in order to discuss 
the numerical impact of the potential problems in practice. We show numbers for both 
the (N)NLLi and (N)NLL2 matching prescriptions, and use E — )■ Note that for 

a fixed soft scale not all the logarithms of /3 are resummed locally in the partonic cross 
section, so the implementations (N)NLLi and (N)NLL2 differ not only by a constant but 
also by logarithmic terms at 0{a'j.). The soft scale is chosen as in the second line of (14.111) . 
For the factorization scale we take /i/ = rrit, for the hard scale fih = 2mt, as discussed 
in Section [2.2.11 The Coulomb scale flc is defined in (12.121) . We estimate the theoretical 
uncertainty as follows: 

Scale uncertainty: We vary all scales Hi (the factorization, soft, hard and Coulomb 
scales) in the interval [/ii/2, 2/ij] around their central values fli. fis and fic are varied 
while keeping the other scales fixed, and fij are allowed to vary simultaneously. 



^For the production of heavier particles, where the threshold region is more important, this ambiguity 
is indeed smaller. For instance, for rrit = 1 TeV and the LHC with y/s = 7 TeV we find p,s — 188 GeV 
expressing the soft function in terms of \og{E / ^g), and fig = 149 GeV using \og{mtf3'^ / iif), so the difference 
is reduced to less than 30%. 



19 



Tevatron 



LHC (x/i = 7 TeV) LHC (x/i = 14 TeV) 



NLL 



6.60 



+0.43+0.00 +0.10+0.50 
-0.65 -0.61 -0.10 -0.44 




692 



+ 169+ 0+28+46 
-180 -85 -28 -41 



NLL2 



6.90 



+0.32+0.00 +0.10+0.52 
-0.41 -0.04 -0.10 -0.47 



157.6 



+23.2+2.3 +4.7+13.5 
-20.2 -0.00 -4.7-12.8 





151.8 



+9.4+5.5 +4.7+13.7 
-4.9-0.0-4.7-13.2 



837 



+69 +36 +28 +57 
-28 - -28 -57 



NNLL2 7.08 



+0.15+0.07 +0.10+0.69 
-0.28 -0.00 -0.10 -0.53 



157.4 



+8.8+5.4+4.7+14.5 
-3.6-0.0-4.7-13.9 



868 



+63 +29 +28 +61 
-21 - -28 -60 



Table 1: Results for Method 1, rrit = 173.3 GeV, all numbers in pb. The four errors refer to 
scale variation, resummation ambiguity, the NNLO constant term, and the RDF+a^ uncertainty. 

imposing the additional constraint 1 < < 4. The errors from varying {fif,fih}, 

Us and fic are obtained by taking the respective maximum/minimum, and are added 
in quadrature. 

Resummation ambiguities: We estimate ambiguities in the resummation procedure by 
taking the difference between the default E — m^/?^ and E = \/§ — 2mt. With the 
second expression for E the soft scales in the first line of ( 14. lip are used. 

NNLO-constant: All constants Cpp) ^ in (14. 5 p are varied by ±\C^J'^^'\. 

PDF+as uncertainty: We estimate the combined error of the NNLL results due to un- 
certainties in the PDFs and the strong coupling using the 90% confidence level eigen- 
vector set of the MSTW08NNLO PDFs [64J and the five sets for variations of 
provided in |69], which corresponds to as{Mz) = 0.1171 ± 0.0034. For the NLL 
prediction the corresponding NLO PDFs are used. 

For the top quark pole mass rrit = 173.3 GeV, we obtain the total top pair production 
cross sections displayed in Table [H where the four errors refer to the four sources of 
uncertainty detailed above. Both, at the NLL and NNLL order, the central values for 
the two matching procedures lie inside the common error band. The total relative theory 
error, given by the sum in quadrature of the scale, resummation and constant uncertainties, 
decreases from NLL to NNLL, as one would expect from including higher orders in the 
logarithmic expansion, and increases from Tevatron to LHC, consistent with the fact that 
the threshold kinematics is parametrically more dominant at Tevatron, and thus matching 
and resummation ambiguities are expected to be smaller. 

The difference of the NLLi results, both to the NNLLi and the complete fixed-order 
NLO predictions, (see Tables IHl and [TO] in Section|5]) are large, especially at the LHC, since in 
NLLi the fixed-order NLO corrections are not included through matching. The difference 
between NLL2 and the fixed order NLO results, as well as the two NNLL implementations, 
is small in comparison. Since part of the difference is due to the use of different PDFs (NLO 
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vs. NNLO), it is instructive to compare predictions using NNLO PDFs also at NLL in order 
to ascertain the genuine effect of tlie liiglier-order corrections. For tlie NLL2 approximation 
we tlien find tlie central values o"^^^^ = 6.67 pb at the Tevatron and cr^^^^ = 148.6 pb 
at the LHC with ^fs = 7 TeV (ag^^^ = 830 pb with ^fs = U TeV). Since in both the 
NLL2 and NNLLi prescriptions the resummed cross section is matched to the fixed-order 
NLO result, their difference shows that the genuine effect of going from NLL to NNLL is 
moderate, and leads to an enhancement of about 3% at the Tevatron and below 1% at 
the LHC with ^/s = 14 TeV. The difference of the two NNLL implementations is of the 
order of 3 — 4%, and is an estimate for the effect of the constant terms at NNLO and the 
difference of including all log /? terms exactly at NNLO or approximately due to the use of 
a fixed /is- 

The total error, obtained by summing in quadrature the four sources of uncertainty, is 
dominated by the scale uncertainty and the error associated with PDF and as variation. 
Despite the ambiguity in the determination of the soft scale (14. lip , the error associated 
to resummation is small at both Tevatron and LHC. The scale uncertainty decreases from 
NLL to NNLL, and is systematically smaller for the second matching option compared to 
the first. 

4.2 Method 2 

For the reasons explained in Section 14. H the use of a fixed soft scale may be problematic 
for top-quark production at the Tevatron or LHC, where the total cross section is not 
dominated by the threshold region, but where it is nevertheless useful to include logarithmic 
higher-order corrections. We therefore propose a procedure that uses a running scale and 
fixes the soft scale only in the low-/3 region. We divide the convolution of the partonic 
cross section with the parton luminosity f l2.16p into two regions using a parameter P^nt, 
chosen such that on the one hand the perturbative expansion in is not spoiled by 
large logarithms in the upper interval /3 > Pent, and on the other hand ambiguities in 
the threshold approximation are small in the lower interval /3 < Pent- The partonic cross 
section is then treated differently in the two regions: 

P < Pcut For Pent chosen sufficiently small (say Pent < 0.5) the threshold expansion 
in P is convergent in the lower interval and the use of the factorization formula (12. ip is 
justified. In this region, the logarithms of P become large and need to be resummed for a 
reliable prediction. The threshold-enhanced contributions are numerically dominant and 
ambiguities in the resummation procedure, like the kinematical ambiguities or the use of 
the different matching prescriptions (14. ip and (14. 2p . should not lead to large numerical 
differences. For the reasons discussed before we use the NNLL2 implementation defined 
in (14. 2p . with E — )■ rritP'^ and a fixed soft scale 

fxf = Km,Pl,. (4.12) 

The choice of the constant is discussed below. 
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/3 > /Scut In the upper interval, the use of the threshold approximation cannot be justified 
a priori. Nevertheless we here adopt the point of view that the inclusion of a subset of 
higher-order corrections is useful even in this region. For Pcut chosen large enough, so that 
the perturbative expansion is not spoiled by large logarithms, the numerical difference of 
the NNLL resummed cross section and the expansion to NNLO and N^LO accuracy is 
expected to be small, so that the result should not depend too strongly on the choice of 
one of these approximations. By default we use again the implementation NNLL2 with 
E — )■ mt/3'^. Since for /9 — >• 1 it is not possible to argue that a fixed soft scale correctly 
includes dominant contributions to the hadronic cross section, we choose a running soft 
scale, 

= ksmt(3\ (4.13) 

in order to include all NNLL contributions to the partonic cross section correctly. In order 
to estimate the ambiguities in this treatment, we compare it to the approximate fixed-order 
cross sections NNLOapp [36], as well as N^LOa and N^LOb defined in Section [3], in each 
case including all lower-order terms up to the indicated order. 

There is a certain tension in the requirements on /3cut- If chosen too large, the cross 
section result becomes too sensitive to the ambiguities in the matching procedure and the 
resulting constant terms at Furthermore, a larger /3cut implies a larger fif, so that 

resummation becomes ineffective. On the other hand, for Pent too small the difference 
among successive orders in the perturbative expansion in the upper interval becomes too 
large (c.f. Figure [1]). This suggests the following procedure to determine /3cnt'- we consider 
the cross section 

au{A^, 5>, /3eut) = OWcut - /?) + - Pent) , (4.14) 

defined using one of the appropriate approximations for the upper and lower intervals, 
A< G {NNLLi,NNLL2} and 5> G {NNLL2, NNLOapp, N^LOa, N^LOb}, and determine 
/Scut such that the difference between the eight different implementations (T(f(yl<, i?>, /3cut) 
becomes minimal. 

We repeat this procedure for values ks = 1,2,4 in the soft scales (14.121) and (I4.13P in 
the two intervals. As default we adopt kg = 2, and obtain the following values for (3cut and 
the corresponding fixed soft scale in the lower interval: 

/3cut(NNLL) : 0.35 (Tevatron), 0.54 (LHC7), 0.55(LHC14), 

/i< = 2mt/32^t : 42 GeV (Tevatron), 101 GeV (LHC7), 105 GeV (LHC14). ^' ' 

It is seen that the /Scut values are in a region where we expect the /3-expansion and the 
perturbative expansion in both to be reasonably reliable, so they satisfy the requirements 
discussed above. We treat the difference to the results for ks = 1,4 as another source of 
the theoretical uncertainty. More details on this procedure and the (3cut values for /c^ = 1, 4 
are given in Appendix [Dl 

To summarize, our default implementation is given by the option NNLL2 defined in (14. 2 p 
using the soft scales (14.120 and (I4.13p . with fc^ = 2 in both intervals and the /Scut-values given 
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Tevatron 



LHC (x/i = 7 TeV) LHC (^/i = 14 TeV) 



NLL2 



7.31 



+0.25+0.30 +0.10+0.57 
-0.03 -0.53 -0.10 -0.54 



172.8 



+14.8+13.0+4.7+15.9 
- 0.8 -14.7-4.7-14.6 





162.6 



+4.2+3.9+4.7+15.4 
-1.9-5.6-4.7-14.7 



896 



+22+18+28 +65 
- 5 -23 -28 -64 



Table 2: Results for Method 2, rrit = 173.3 GeV, all numbers in pb. The four errors refer to 
scale variation, resummation ambiguity, the NNLO constant term, and the RDF+a^ uncertainty. 

in (14.151) . For the remaining scales we again take the default values flf = rrit, jlh = 2mj 
and the Coulomb scale as defined in (I2.12p . We estimate the remaining uncertainties as 
following: 

Scale uncertainty: as for Method 1, excluding the soft scale. The latter is effectively 
replaced by the variation of ks. 

Resummation ambiguities: we consider three different sources for ambiguities: i) The 
difference between the default setting E = rritfi'^ compared to E = a/S — 2mt in the 
NNLL2 implementation, ii) the difference between the NNLL2 implementation for 
the soft scale choices kg = 1,4 (and the corresponding /3cut values) to the default 
choice ks = 2, iii) the envelope of the 8-different approximations (I4.14p for variations 
of /3cut by 20% around the default values (I4.15P for k^ = 2. The resulting errors are 
added in quadrature. 

NNLO-constant, PDF+a^ uncertainty: estimated as in Method 1. 

Our results for Method 2 are shown in Table El where the four errors refer to the four 
sources of uncertainty as detailed above. The separate sources of resummation ambiguity 
are shown in more detail in Table [3l from which one can see that the variation of kg and 
/Scut (i-e. the soft-scale variation) dominate the error, while the ambiguity due to the use of 
E = rritP'^ is negligible. To see the genuine effect of NNLL resummation it is again useful 
to consider the NLL2 results using the NNLO PDFs for which we obtain aJJ^^^ = 7.02 pb 
at the Tevatron and aJJ^^^ = 162.2 pb at the LHC with y/s = 7 TeV (aJJ^^^ = 901 pb at 
the LHC with y/s = 14 TeV). Therefore the effect of the NNLL corrections on the central 
values is small, about 3% at the Tevatron and in the per-mille range at the LHC, showing 
a good convergence of successive orders of the logarithmic approximations. 

The NNLL2 results in Table |2] are about 3% larger than those obtained in Table [1] using 
Method 1, while for the NLL2 results the differences are of the order of 10%. However the 
results from the two methods agree within the estimated resummation ambiguity. Since 
the soft scales corresponding to (3cnt in ( 14.15P lie in the range of the fixed soft scales used 
for Method 1 (14. lip , the differences between the two methods are smaller than might have 
been expected, and serve as an estimate of the effect of treating logarithmic corrections 
exactly for large /3, as in Method 2, or approximately through a fixed soft scale, as in 
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Collider i) E ii) ks iii) Pent 



Tevatron 




+0.01 
-0.00 


r+0.l%x 

L-0.0%J 


+0.10 
-0.10 


r+1.4%1 
1-1.4%/ 


+0.17 
-0.18 


r+2.4%x 
1-2.5%/ 


LHC (v^i = 


7 TeV) 


+1.0 
-0.0 


r+0.6%1 
L-0.0%J 


+0.8 
-2.2 


r+0.5%T 
1-1.4%/ 


+3.7 
-5.1 


r+2.3%1 
1-3.1%/ 


LHC (v^ = 


14 TeV) 


+5 
-0 


r+0.6%x 
L-0.0%J 


+ 

-10 


r+0.0%x 

1-1.1%/ 


+17 
-21 


r+1.9%X 
1-2.3%/ 



Table 3: Different resummation ambiguities of the NNLL2 cross section for rrit = 173.3 GeV. 
Absolute errors are in pb. Relative errors are given in brackets. 



Method 1. Since the results of both (very different) methods for the treatment of the 
soft scale agree within errors, we are confident that we are realistically estimating the 
resummation ambiguities. 

The total relative theoretical error (not including the PDF+as uncertainty) is shown 
in Table S] for the NLL2 and NNLL2 cross sections, and for the NNLOapp result (see 
Section O Tables |H] and [TDl for the absolute values). Both for the Tevatron and the 
LHC, the theoretical error decreases when higher-order logarithmic corrections are added 
to NLL2. However at the Tevatron the smallest uncertainty is obtained for NNLOapp, 
while NNLL2 yields the smallest error estimate for LHC. To understand this somewhat 
surprising behaviour at the Tevatron, in Table O we study separately the results for the 
quark- ant iquark and gluon-gluon channels. It is seen that the increase of the error estimate 
for the quark-antiquark induced channel from NNLOapp to NNLL2 occurs both at the 
Tevatron and the LHC. On the contrary, at the LHC the error of the dominant gluon 
channel is reduced at NNLL2, which leads to the observed reduction of the total theoretical 
uncertainty. At the Tevatron a reduction of the error in the gluon channel is observed as 
well, if the same value /3cut = 0.54 as for LHC is used. It is interesting to recall that 
a comparison of the NLO singular terms with the exact NLO results (see e.g. [381139]) 
shows that, while the exact partonic cross section for the gg channel is approximated 
reasonably well by its threshold approximation over the whole /3 range, for the qq channel 
the approximation breaks down at /3 ~ 0.3. This points to a more problematic behaviour 
of the quark-antiquark channel. Therefore one might conclude that the small error of the 
NNLOapp prediction at the Tevatron is an underestimate of the true uncertainty, due to 
the particular behaviour of the dominant quark-antiquark channel]^ 

To summarize, by using a running soft scale within Method 2, we implement the princi- 
ple of scale separation in effective theories at the level of the partonic cross section, rather 

^This observation could suggest a separate /3cut choice for both channels, however due to the small 
effect on the central value we continue to use the same values. 

^Another indication of this comes from a study of the choice fir — 2mt for the central value of the 
renormalization scale in NNLOapp, that is identical to the central value of /i^ in the NNLL results. For 
this choice, the scale-uncertainty increases moderately to 2% for NNLOapp at the LHC with -y/s = 7 TeV, 
while the error at the Tevatron increases to ^g g% and becomes larger than for NNLL2. 
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Tevatron 



LHC ( = 7 TeV) LHC (^/i = 14 TeV) 



NLL2 {+5.5%, -7.4%} {+11.7%, -8.9%} {+11.6%, -8.0%} 

NNLOapp {+3.8%, -4.9%} {+7.7%, -7.4%} {+8.6%, -7.7%} 
NNLL2 {+4.2%, -6.5%} {+4.6%, -4.6%} {+4.5%, -4.1%} 

Table 4: Relative error of NLL2, NNLOapp and NNLL2 at the LHC. Scale uncertainty, resum- 
mation ambiguity and the error from the NNLO constant are added in quadrature. 

Channel Tevatron LHC {^/s = 7 TeV) 

NNT n nn 10 +0.21 +0.07 /+3.6%\ 90 7+2.6+0.4 / +9.2% \ 

iNiNLUapp qq O-J-^ -0.27 -0.07 I -4.5%; ' -2.9 -0.4 | -10.2% J 

ATATT T „77 01 +0.23+0.12+0.07 /+4.3%\ 90 7 +2.4 +0.3 +0.4 / +8.5% \ 

iMMjl^2 qq O.di_o.4i_o.20-0.07 |-7.3%J ' -3.0 -1.9 -0.4 | -12.5% J 

NNTO nn n Qf; +0-10 +0-02 / +10.7% \ 1 00 9 +17.4 +4.3 / +13.5% \ 

IMiMljWapp yy "J-y^ -0.09 -0.02 I -9.7% J -'-'-''-'•^-13.8-4.3 |_-10.9%J 

ATATTT n Qo +0.02 +0.12 +0.02 /+12.6%\ 10/17+8.2 +3.6+4.3 / +7.4% \ 

iNiNijij2 yy U.y«_o.o8 -0.05 -0.02 \ -9.S% / ' -12.5 -3.9 -4.3 \-10.2%J 



ATATTT „„ n QQ +0-06 +0.05 +0.02 1+8.1% 1 10/17+8.2 +3.6+4.3 1+7.4% 

iNiNijij2 yy '-'•yy -0.08 -0.06 -0.02 1 -10.3%! -^'-'^■'-12.5-3.9-4.3 ) -10.2% 

(/3cut = 0.54) ^ ^ ^ ' 



Table 5: Results for the qq- and yy-channe\, Method 2, rrit = 173.3 GeV, all numbers in pb. 
The errors refer to scale variation, resummation ambiguity (for NNLL2 only) and the NNLO 
constant term. The numbers in brackets denote the total relative error. Unless stated otherwise, 
the default values of /3cut, Eq. (14.150 . are used. 



than at the level of the hadronic cross section as in Method 1. Threshold logarithms are 
resummed locally rather than in an average sense in the region above Pcut, while problems 
related to very small soft scales for small /3 are avoided. The sensitivity to the precise value 
of /3cut is moderate and included in our error estimate. In Section [5] we use this method as 
our default implementation. 



4.3 Effect of Coulomb resummation 

In this section we would like to comment on the size of the bound-state contributions and 
Coulomb resummation, which were not (or only partially) included in previous works on 
tt threshold resummation. 
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Toponium-like bound states may form below the nominal production threshold due 
to the strong Coulomb attraction in the colour- singlet channel. In practice, the bound- 
states are smeared into a broad resonance due to the rapid t — )■ bW decay. For an inclusive 
observable such as the total cross section the resonance effect can be computed by neglecting 
the top decay width and summing instead over the series of would-be toponium bound- 
state poles in the Coulomb Green function. The technical aspects are discussed in detail 
in Appendix |Bl In Table |6] we show the corresponding contribution to the cross section, 
denoted by BS, which are always included in the NNLL results we present. The first error 
denotes the sum in quadrature of scale and resummation uncertainties, determined as from 
Method 2, and the second one the PDF+a^ error. The bound-state contributions, which 
contribute first at 0{al), are rather small and make up less than 0.5% of the total cross 
section. The large theoretical error might be due to the more singular behaviour of the 
convolution of the bound-state correction with the soft function and the PDFs (c.f. (1B.6P ). 
which makes the cross section more sensitive to variations of /i/ and kg. 

Next, we discuss the effects due to the resummation of Coulomb corrections in the 
continuum [E > 0). To switch off the Coulomb effects in the resummation formula, we 
set the Coulomb function (12. lip in Eq. 12.11 to its first-order term in the expansion in ag 
(see (EH)), 



We then consider the quantities 



5Cbi = a^NLL, _ NNLLi 



tnv 



5Cb2 = a^NLL. _ ^NNLL, ^ ^^^^^^ 

where (j'^^'^'^* denotes the total resummed cross section with the full Coulomb function (but 
excluding bound-state effects), matched according to prescription (14. ip or (14. 2p . and cr^iv'"'"* 
the analogous quantity for the trivial potential function (I4.16p . In this case, the expansion 
of 

Ctriv^^ to NLO or NNLO order is subtracted in the matching prescriptions (14. ip or (14. 2p . 
As a consequence of the matching to the fixed-order results, 5Cbi contains only terms of 
order or higher. In particular, it contains all the purely Coulomb contributions, and 
the interference of Coulomb corrections with hard and soft terms, at NNLO. 5Cb2 contains 
only terms of order or higher, and gives a measure of the effect of Coulomb resummation 
beyond NNLO. 

In Table [6] we show the central values for 5Chi and 5Cb2- We do not show errors, which 
are large, due to 5Cbi^2 being defined as a difference. From the numbers in the table, it 
is clear that the dominant effect of Coulomb resummation is accounted for by the terms 
already included in NNLOapp, which coincide with the difference (5Cbi — (5Cb2. They give a 
non-negligible correction at both Tevatron (~ 2%) and LHC (~ 1 — 2%). Beyond NNLO, 
Coulomb corrections, and their interference with soft and hard resummation, is very small 
(^ 0.5%), as can be inferred from the numbers for 5Cb2. 
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Tevatron 



LHC (v^ = 7 TeV) LHC (v^ = 14 TeV) 



BS 





+ 1.7+0.5 
-1.6-0.4 



6Ch 



1 



0.140 



2.81 



7.8 



(5Cb2 



-0.052 



0.13 



-0.3 



Table 6: Bound-state and Coulomb contributions for rrit = 173.3 GeV, Method 2. All numbers 
in pb. 

4.4 Comparison to other NNLOapp and NNLL predictions 

Previous results for the total top cross section beyond NLO and NLL |^[66] use theoretical 
input equivalent to the NNLOapp result discussed in the present paper (and presented in de- 
tail in the subsequent section), but the complete NNLL threshold resummation (including 
Coulomb effects) performed here is new. 

There exists a complementary approach to resummation, which is based on resumming 
logarithms in various kinematical limits different from the partonic threshold /3 — )■ such 
as pair invariant mass (PIM) kinematics [37|l38] to compute the invariant mass distribution, 
and one-particle inclusive (IPI) kinematics [1T|H2] for rapidity and transverse- momentum 
distributions. Results for the total cross section are then obtained by integrating PIM or 
IPI differential cross sections, and do not include the higher-order Coulomb corrections, 
which, as discussed in section H73| amounts to an effect of —0.13 pb at the Tevatron and 
+3.5 pb (+11 pb) at the LHC with = 7 TeV (v^ = 14 TeV). On the other hand, these 
calculations keep certain sets of power-suppressed terms in /3 that can be non-negligible, 
since the total cross section is dominated by contributions with /3 > 0.3. Therefore these 
approaches complement our treatment using the partonic threshold limit, and differences 
in the predictions can indicate systematical uncertainties of resummation. In Table [7] we 
compare the results from these approaches, which are either based on NNLOapp expansions 
or on NNLL resummations to ours. 

The approximations based on IPI and PIM kinematics are themselves subject to large 
ambiguities, analogous to the ambiguity between E = y/I — 2mt and m^/?^ in the threshold 
limit. In [381 |l2] if is argued by comparison to the known fixed-order NLO results that 
a particular form of the singular distributions arising in the SCET-based resummation 
formalism is preferred over those used conventionally (e.g. [T7jl^ ). This ambiguity between 
the conventional IPI approximation and the IPIscet implementation of [42j is reflected 
in the different predictions in the first two rows of Table [3 It is interesting to note that 
these kinematical ambiguities appear to be larger in IPI and PIM kinematics than in the 
threshold approximation where we only found a small ambiguity (c.f. Table |3]). 

From the results in Table [7] it is seen that our results agree within the quoted errors 
with those of Kidonakis jH] in IPI kinematics. While the tension with the predictions of 
Ahrens et. al. [381H2] is larger, the results for the LHC agree within the given errors, if the 
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Tevatron LHC (v^ = 7TeV) LHC (v^ = 14TeV) 



± > ± > ij w app 
app 

NNLO^j^^scET 


(Ref. m) 
(Ref. m) 
(Ref. [38]) 


7 08+0-20 
' •U«_o.24 

fi fiq+000 

U.UO_Q 27 
U.UZ_Q 4Q 


1631^ 
1551^ 

1551^ 


9201^° 

85ltf 
860;^^ 


NNLL^^^scET 


(Ref. m) 
(Ref. [38]) 


fi c^c;+o-i6 

"•^'-'-0.14 

fi 46+0-18 

"•^"-0.19 




8241^1 
81111^ 


NNLL2 


this work 


7 99+0.31 
' -^^-0.47 


163tl 


896;^° 



Table 7: Comparison to previous NNLOapp and NNLL predictions for the central value fif = rrit 
and MSTW08NNLO PDFs, all numbers in pb. In the results of [4T] mt = 173 GeV, in those 
of [3B1|42] rrit = 173.1GeV. Only theory errors, excluding PDF and errors, are shown. The 
NNLL2 theory errors from table [2] have been added in quadrature. 

size of the Coulomb corrections included in our results is taken into account. The size of 
the theory uncertainties quoted by the different groups are comparable at the LHC (except 
for the very small error assig ned to NNLO!,ppSCET)^ however our theory error includes an 
estimate of the missing NNLO constant and the resummation ambiguities, which are not 
included by the other groups. At the Tevatron, the agreement with the results of [38lll2] is 
worse, in particular for the resummed results. In this case our estimate of the theory error 
is more conservative than that of the other groups^ The larger discrepancies between 
results obtained using different kinematical approximations might again point to a worse 
behaviour of the threshold approximation for the quark-antiquark channel dominant at the 
Tevatron, as discussed already at the end of Section 14.21 

5 Detailed numerical results 

In this section we present numerical results for the combined NNLL resummation of soft 
and Coulomb effects for the total ti cross section at hadron colliders. This requires a choice 
of the parton distribution functions used in the convolution of the partonic cross section 
with the parton luminosity. Available fits for PDFs at NNLO accuracy are the sets of 

""^"^The quoted results from |41j for the Tevatron include an independent variation of /i/ and fir while 
the NNLOapp results from [38l|42| set /if — Hr- Adding separate variations of both scales in quadrature 
and averaging the IPIscet and PIMgcET kinematics, Ahrens et.al. [70] recently obtained a slightly in- 
creased error estimated a^^ = 6.63+q'4^ pb. It is interesting to note that the difference between the 
two kinematics is larger for the resummed cross sections than for the NNLOapp-results combined in the 
prediction from |70| . 
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MSTW2008NNLO P], JR09 [TT], ABKM09 [72], HERAPDFl.O [73] and NNPDF2.1 [71]. 
Recent discussions of the impact of the differences among the different PDF fits on the 
top-pair cross section can be found in [74l[75], where it was shown that the MSTW and 
NNPDF sets yield consistent results, whereas there are larger differences between the other 
sets, partly due to different values of the strong coupling constant obtained by the different 
groups. For our results we use the MSTW2008NNLO PDF sets at 90% confidence level 
and the associated value of the strong coupling constant |64j . 

«.(M|) = O.lUltro'olt ■ (5.1) 

Besides numbers for the NNLL cross section, we also present the NLO result [TT1I65] and 
NNLOapp, which adds to NLO the threshold expansior0 of the NNLO correction [36] 
(see Appendix [C] for a summary of the analytical expressions) to illustrate the signifi- 
cance of higher-order results. For NLO we use the NLO PDFs and the associated value 
as(M|) = 0.1202lQ;QQ3g. In all three cases, we determine the joint PDF+as uncertainty 
using the method described in [69] . Following [16] the theoretical uncertainty of the fixed- 
order NLO and NNLOapp approximations is determined by setting the central values of the 
renormalization and factorization scales, {fir, /i/}, to rrit and varying both scales simulta- 
neously in the interval [/ii/2, 2/ij], imposing the additional constraint 1/2 < fj^r/f^f < 2. 

For the NNLL resummed cross section we choose the /3-dependent soft scale (14.12^ . 
f l4.13p . and the matching condition (14. 2p . denoting the corresponding result by NNLL2. The 
error estimate follows the procedure detailed in Section ^^21 The errors from scale variation 
and the resummation ambiguities are added in quadrature, while the uncertainty from 
variation of the unknown NNLO constant is given as a separate error both for NNLL2 and 
NNLOapp. The default value of the pole mass of the top is chosen as rrit = 173.3 GeV [76] . 
Results for rrit = 165 . . . 180 GeV are given in Table |9] for Tevatron, and in Tables [TT] and 
[T2] for LHC In Sections 15.11 and 15.21 we also compare our results to recent measurements 
at Tevatron and LHC, while in Section 15.31 we illustrate the impact on the extraction of 
the top mass from measurements of the cross section. 

5.1 Tevatron 

The cross section results for Tevatron kinematics are summarized in Table El The differ- 
ence between NNLL2 and NLO amounts to about 8% of the NLO result. This correction 
is an interplay of genuine resummation effects, which amount to about 12% of the NLO 
result, and of switching from NLO to NNLO PDFs, which lowers the NLO cross section 
from the value given in the table to 6.46 pb. As evident from Table [H the bulk of the 
resummation corrections are accounted for by the fixed-order soft and Coulomb terms at 
0{ai) (NNLO app), which give a 9.5% correction to the NLO result, though higher-order 

^^Note that we apply a strict threshold approximation for NNLOapp, i.e. we discard all terms scaling 
as 0{a1/3'^) or higher with respect to the leading Born cross section. Consequently, our central value for 
NNLOapp differs, albeit slightly, from the one obtained e.g. using HATHOR, ^66J, since there ©(a^/?^) 
terms from the gq-cha,imel are included. Furthermore, terms of the form const, x In (/i//mt), although 
known exactly at NNLO, are not included, in the spirit of a strict threshold expansion. 
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NLO 



NNLL2 



+0.36+0.51 7 nfi +0.25+0.10+0.69 7 99 +0.29+0.10+0.71 

-0.75-0.45 '•'-'"-0.33-0.10-0.53 '-^^-0.46-0.10-0.55 



Table 8: The total top-pair cross section (in pb) for nit = 173.3 GeV at the Tevatron. The 
first set of errors refers to scale variation (scale variation+resummation ambiguities for NNLL2), 
the last to PDF+«s error. The second set of errors for NNL0app/NNLL2 arises from variations 
of the unknown NNLO constant term. 



o-„[pb] 



Tevatron: Theory Uncertainty 




r m,[GeV] 

166 168 170 172 174 176 178 180 

o-,r[pb] 



Tevatron: Theory Uncertainty 




m,[GeV] 

166 168 170 172 174 176 178 180 



Figure 3: Dependence of the total cross section on the top mass at the Tevatron. The bands 
denote the total theory uncertainty, i.e. the sum (in quadrature) of the scale and resummation 
ambiguities, and the estimate of the NNLO constant but not the PDF + error. 
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Table 9: Total cross sections in pb at the Tevatron for rrit = 165 . . . 180 GeV. The errors 
denote the scale variation (scale variation+resummation ambiguities for NNLL2), the NNLO 
constant variation (for NNLOapp and NNLL2) and the PDF+a^ error. 



rrit [vjeVJ 


ATT r\ 




NNLOapp 


ATATT T 


165 


8.70 


+0.48+0.66 
-0.98-0.61 


9.17 


+0.33+0.13+0.93 
—0.44—0.13—0.71 


Q 00 +0.39+0.13+0.96 
y'-JO -0.60-0.13-0.73 


166 


8.42 


+0.46+0.64 
-0.95-0.58 


8.88 


+0.32+0.13+0.89 
—0.42—0.13—0.68 


Q (-,0+0.37+0.13 + 0.92 

y-Uo -0.58-0.13-0.71 


167 


8.16 


+0.45+0.62 
-0.92-0.56 


8.60 


+0.31+0.12+0.86 
—0.41—0.12—0.66 


„ op, +0.36+0.12+0.89 
°-°U_o.56-0.12-0.68 


168 


7.90 


+0.43+0.60 
-0.89-0.54 


8.33 


+0.30+0.12+0.83 
— 0.39—0.12 — 0.63 


Q rr) +0.35+0.12+0.86 
o-^-^ -0.54-0.12-0.66 


169 


7.65 


+0.42+0.58 
-0.86-0.52 


8.07 


+0.29+0.11+0.80 
—0.38—0.11—0.61 


Q r)(:> +0.34+0.11+0.83 
°-^t)_o.53_o.ii_o.64 


170 


7.41 


+0.41+0.56 
-0.84-0.51 


7.82 


+0.28+0.11+0.78 
— 0.37—0.11 — 0.59 


riri +0.32+0.11+0.80 
0-UU_o.5i_o.ll-0.61 


171 


7.18 


+0.39+0.54 
-0.81-0.49 


7 58 


+0.27+0.11+0.75 
-0.36-0.11-0.57 


w +0 31+0 11+0 77 
' ■ -0.49-0.11-0.59 


172 


6.96 


-|-U.oo-rU.Oo 

-0.78-0.47 


7.35 


+0.26+0.10+0.72 
-0.35-0.10-0.55 


' •"-'^-0.48-0.10-0.57 


173 


6.74 


+0.37+0.51 
-0.76-0.46 


7.12 


+0.25+0.10+0.70 
-0.34-0.10-0.53 


7 9Q +0.29+0.10+0.72 
' -0.46-0.10-0.55 


174 


6.54 


+0.36+0.50 
-0.74-0.44 


6.91 


+0.24+0.09+0.67 
-0.32-0.09-0.51 


7 r, 7 +0.28+0.09+0.70 
' ' -0.45-0.09-0.53 


175 


6.34 


+0.34+0.48 
-0.71-0.42 


6.70 


+0.23+0.09+0.65 
-0.31-0.09-0.49 


c ocr +0.27+0.09+0.67 
U.OO_Q 44_Q Q9_Q 52 


176 


6.14 


+0.33+0.47 
-0.69-0.41 


6.49 


+0.22+0.09+0.63 
-0.30-0.09-0.48 


c f^A +0.26+0.09+0.65 
"•"^-0.42-0.09-0.50 


177 


5.96 


+0.32+0.45 
-0.67-0.40 


6.30 


+0.22+0.09+0.61 
-0.30-0.09-0.46 


f^AA +0.26+0.09+0.63 
"•^^-0.41-0.09-0.48 


178 


5.78 


+0.31+0.44 
-0.65-0.38 


6.11 


+0.21+0.08+0.59 
-0.29-0.08-0.45 


f- 9 c- +0.25+0.08+0.61 
"•^'-'-0.40-0.08-0.46 


179 


5.60 


+0.30+0.43 
-0.63-0.37 


5.92 


+0.20+0.08+0.57 
-0.28-0.08-0.43 


nfi -I-0-24+0. 08+0.59 
"•""-0.39-0.08-0.45 


180 


5.43 


+0.29+0.42 
-0.61-0.36 


5.75 


+0.19+0.08+0.55 
-0.27-0.08-0.42 


pr OQ +0.23+0.08+0.57 
^•°" -0.37-0.08-0.43 



contributions (from NNLOapp to NNLL2) amount to a non-negligible 3.5%. The total the- 
oretical error of the cross section is reduced at both NNLOapp and NNLL2 compared to the 
NLO result, and amounts to {+4.2%, -6.5%} for the NNLL2 resuh, and {+3.8%, -4.9%} 
for NNLOapp. As discussed in Section l42l the counter-intuitive increase of the uncertainty 
of the resummed result compared to the fixed-order NNLO is related to the behaviour 
of the quark-antiquark induced partonic subprocess. Interestingly, the relative PDF+a.; 
uncertainty increases when going from NLO to NNL0app/NNLL2, an effect which is also 
observed at the LHC and has been found in [75] as well. The ambiguity arising from 
the unknown constant (second error) is small compared to the scale, resummation and 
PDF+as uncertainties, which dominate the total error at NNLL2. 

The dependence of the total cross section on rrit is plotted in Figure [3] and given explic- 
itly in Table [9l The relative theoretical errors in the mass range rrit = 165 . . . 180 GeV differ 
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v/i NLO 



NNLL2 



7 



14 




891 




+71+28+64 
-63-28-63 




162.6 



+5.7+4.7+15.4 
-5.9-4.7-14.7 



Table 10: Predictions for the total cross section (in pb) for rrit = 173.3 GeV at the LHC. The 
errors are defined as in Table [HI 

only on the permille level from those for rrit = 173.3 GeV, while the PDr+a^ uncertainty 
decreases slightly for larger masses, resulting in a practically constant overall uncertainty 
for NNLL2 of tllii at TJit = 165 GeV and l^^^g'" at = 180 GeV. 

The NNLL2 results are in good agreement with measurements performed at the Teva- 
tron. The DO experiment obtains au = 7.56^Q5g pb from combining measurements of the 
dilepton and lepton plus jets final states with up to 5.6fb~^ of data [1]. The CDF collabo- 
ration quotes au = 7.50^Q;4g pb from a combination [2] of the dilepton, the lepton plus jets 
and the all-hadronic channel using up to 4.6fb~^ of data. A value of au = 7.70to;52 pb has 
been obtained from a measurement of the ratio of the top-pair and Z-boson production 
cross sections [1]. In these measurements a top mass rrit = 172.5 GeV has been used in the 
kinematical reconstruction of the top-events. 



Our results for the LHC with centre-of-mass energies y/s = 7 and 14 TeV are shown in 
Table [TOl In Figure S] we plot the dependence on ^/s of the NNLL2 cross section. The 
difference between the NNLL2 and NLO results is smaller than at Tevatron, and amounts 
to 3% (1%) at ^/s = 7 TeV {^/s = 14 TeV). Again, the genuine effect of resummation 
is partly hidden by the switch to NNLO PDFs for the resummed results. Using NNLO 
PDFs for the NLO result would yield 148.9 pb and 837 pb for the NLO cross section at 
7 and 14 TeV, which differ from the NNLL2 results by 9% and 7%, respectively. As for 
the Tevatron, the bulk of the corrections beyond NLO come from the O(a^) terms, with 
higher order contributions at the 1% level. 

The total theoretical error is largely reduced from NLO to NNLOapp, and from NNLOapp 
to NNLL2 as seen from Table [TOl For NNLL2 it amounts to a {-|-4.5%, —4.6%} residual 
uncertainty at y/s = 7 TeV, and {-|-4.5%, —4.1%} at i/i = 14 TeV. Contrary to the Teva- 
tron, the uncertainty arising from the variation of the unknown constant is comparable to 
the remaining theoretical error. This is a consequence of our estimate \C^J>^^' \ = {C^J r)"^^ 
see Eq. fl4.9p . which gives larger values at LHC, where the gg channel is dominant. At 
NNLL2 the dominant source of uncertainty is the PDF-t-a^ error, which is ±9%. The total 
relative uncertainty, given by the sum in quadrature of all sources of error, is reduced from 
± 15.5% (NLO) to approximately ± 10.5% (NNLL2) at 7 TeV . A similar picture holds 
for 14 TeV, with a smaller total relative error, ± 8.5%. 



5.2 LHC 
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Figure 4: Dependence of the total cross section on y/s at the LHC. The central line is obtained 
with our NNLL2 prediction, while the band corresponds to the sum (in quadrature) of theory 
and PDF+a, errors. 
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Figure 5: Dependence of the total cross section on the mass at the LHC with y/s = 7 TeV. 
The bands denote the total theory uncertainty, i.e. the sum (in quadrature) of the scale and 
resummation ambiguities, and the estimate of the NNLO constant but not the PDF + error. 
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Table 11: Total cross sections in pb at the LHC (^i = 7 TeV) for nit = 165 . . . 180 GeV. 
The errors denote the scale variation (scale variation+resummation ambiguities for NNLL2), 
the NNLO constant variation (for NNLOapp and NNLL2) and the PDF+a^ error. 



mt [GeV] 


NLO 




NNLO 


app 


NNLL2 


165 


203.9 


-|-Z0.0~T"1 1 .0 

-27.4-16.7 


207.6 


-f-iO.UT-D.i-riy.o 
-14.3-6.1-18.6 


209.5 


+7.3+6.1+19.6 
-7.5-6.1-18.8 


166 


197.7 


-1-94 74-1 7 
-26.5-16.2 


201.2 


4-1 A c;-Lf; Q_L1 Q c 

-13.9-5.9-18.0 


203.1 


+7.1+5.9+19.0 
-7.3-5.9-18.2 


167 


191.6 


-LOQ Q_l_l S 

-25.7-15.7 


195.1 


4-14 Ci-L^ 7-l-1fi 9 

-13.4-5.7-17.5 


196.9 


+6.9+5.7+18.5 
-7.1-5.7-17.7 


168 


185.8 


-LO*^ 1 _L1 q 

^^0. -LU. 
-24.9-15.3 


189.2 


4-1 '^4-'^ ^-\-^ 7 7 

-13.0-5.5-17.0 


190.9 


+6.7+5.5+17.9 
-6.9-5.5-17.2 


169 


180.2 


-1-99 4-1-1 ^ A 

-24.2-14.9 


183.5 


4_1 Q 1 4-f; 44_1 7 9 

-12.6-5.4-16.5 


1 8^ 2 


+6.5+5.4+17.4 
-6.7-5.4-16.7 


170 


174.7 


_l91 7_L1 Q 
-f-z i . 1 -f-io.o 

-23.4-14.4 


177.9 


4_19 74_f; 94_1fi 7 
-f-iZ. i-t-O.Z-f-lD. ( 

-12.2-5.2-16.0 


1 79 6 


+6.3+5.2+16.9 
-6.5-5.2-16.2 


171 


169.5 


-22.7-14.0 


172.6 


4-19 '^4_f; flJ-lfi 9 

-f- i z . -f- . u -r ± D . z 
-11.8-5.0-15.6 


17/19 


+6.1+5.0+16.5 
-6.3-5.0-15.8 


172 


164.4 


+20.3+14.5 
-22.0-13.6 


167.5 


+11.9+4.9+15.8 
-11.4-4.9-15.1 


169.0 


+6.0+4.9+16.0 
-6.1-4.9-15.3 


173 


159.6 


+19.7+14.0 
-21.4-13.3 


162.5 


+11.5+4.7+15.4 
-11.0-4.7-14.7 


164.0 


+5.8+4.7+15.6 
-6.0-4.7-14.9 


174 


154.8 


+19.1+13.6 
-20.7-12.9 


157.7 


+11.1+4.5+14.9 
-10.7-4.5-14.3 


159.2 


+5.6+4.5+15.1 
-5.8-4.5-14.4 


175 


150.3 


+18.5+13.3 
-20.1-12.5 


153.1 


+10.8+4.4+14.5 
-10.4-4.4-13.9 


154.5 


+5.5+4.4+14.7 
-5.6-4.4-14.0 


176 


145.9 


+17.9+12.9 
-19.5-12.2 


148.6 


+10.4+4.3+14.1 
-10.0-4.3-13.5 


150.0 


+5.3+4.3+14.3 
-5.5-4.3-13.6 


177 


141.6 


+17.4+12.5 
-19.0-11.9 


144.3 


+10.1+4.1+13.7 
-9.7-4.1-13.1 


145.6 


+5.2+4.1+13.9 
-5.3-4.1-13.3 


178 


137.5 


+16.9+12.2 
-18.4-11.5 


140.1 


+9.8+4.0+13.4 
-9.4-4.0-12.7 


141.4 


+5.0+4.0+13.5 
-5.2-4.0-12.9 


179 


133.5 


+16.3+11.8 
-17.9-11.2 


136.1 


+9.5+3.9+13.0 
-9.2-3.9-12.4 


137.4 


+4.9+3.9+13.2 
-5.0-3.9-12.5 


180 


129.7 


+15.9+11.5 
-17.3-10.9 


132.2 


+9.2+3.8+12.6 
-8.9-3.8-12.0 


133.4 


+4.7+3.8+12.8 
-4.9-3.8-12.2 



The dependence of the cross section on for y/s = 7 TeV in different approxima- 
tions is plotted in Figure |5l Additionally, numerical results for nit = 165 . . . 180 GeV for 
NLO, NNLOapp and NNLL2 at the LHC can be found in Tables [ID and [121 As for the 
Tevatron, the relative theory error is constant to a good accuracy in the considered mass 
interval, while the PDF-t-a^ error increases slightly, resulting in a practically constant over- 
all uncertainty ranging from to tloU at = 7 TeV and from to at 

= 14 TeV. 

Recently the LHC collaborations presented measurements of the total tt cross section 
that approach the experimental accuracy of the Tevatron. Utilizing kinematic information 
of lepton plus jets events, the ATLAS experiment obtains [5] att = 179.0±9.8(stat-|-syst) ± 
6.6(lumi) = 179.0 ± 11.8 pb using up to 0.7fb~^ of data, while the CMS collaboration 
obtains [6] att = 154 ± 17(stat-|-syst) ± 6(lumi) pb from a combination of lepton plus jet 
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Table 12: Total cross sections in pb at the LHC {^/s = 14 TeV) for nit = 165 .. . 180 GeV. 
The errors denote the scale variation (scale variation+resummation ambiguities for NNLL2), 
the NNLO constant variation (for NNLOapp and NNLL2) and the PDF+as error. 



rrit [GeV] 


NLO 




NNLOapp 


NNLL2 


165 


1111 


-4-1 '?fi_Lfi9 

-133-73 


-1 -, -, r-, 4-Q1 4_'^fi4_Sn 
1 1 1 Q -ryi-rOO-|-c5U 

-L-L-L^ -80-36-78 


1124 


+36+36+80 
-30-36-79 


166 


1080 


-129-71 


,^,^^-i-SR4_'^'^4-77 


1093 


+35+35+78 
-29-35-77 


167 


1051 


_1_1 9e_i_7Q 

-126-69 


- ^ p- w 4_Sfl4_'i/l4_7c; 

lUO ( _76_34_74 


1064 


+34+34+76 
-28-34-75 


168 


1022 


+124+76 
-122-67 


1 r^^r^ 4_S'H4-'^'^4^7'^ 

±UZJ_73_33_72 


1035 


+33+33+74 
-27-33-73 


169 


994 


4-1 91 -1-7'^ 
-f- 1 

-119-65 


^ 4_H1 4_'?94_79 

lUUl _7i_32_7i 


1007 


+32+32+72 
-27-32-71 


170 


967 


4-1 1 74-71 

-116-63 


, 4_7Q4_'i1 4_7n 

^'^-69-31-69 


979 


+31+31+70 
-26-31-69 


171 


941 


IT 1 /I 4-70 

-113-62 


-67-30-67 


yoo 


+31+30+69 
-25-30-68 


172 


916 


+111+68 
-110-60 


Q99 +74+29+66 
^^^-65-29-65 


928 


+30+29+66 
-25-29-66 


173 


892 


+108+66 
-107-58 


OQO +72+28+64 
o-'o -64-28-64 


903 


+29+28+65 
-24-28-64 


174 


868 


+105+64 
-104-57 


07/I +70+28+63 
0' ^-62-28-62 


879 


+28+28+63 
-23-28-63 


175 


845 


+102+63 
-101-55 


oci +68+27+61 
-60-27-60 


856 


+27+27+62 
-23-27-61 


176 


823 


+99+61 
-98-54 


000 +66+26+60 
-58-26-59 


834 


+27+26+60 
-22-26-59 


177 


801 


+96+59 
-96-53 


007+64+25+58 
°^ ' -57-25-57 


812 


+26+25+59 
-22-25-58 


178 


780 


+94+58 
-93-51 


78fi +62+25+57 
'00 -55-25-56 


791 


+25+25+57 
-21-25-57 


179 


760 


+91+56 
-91-50 


ycc +61+24+55 
' "J "J -54-24-55 


770 


+25+24+56 
-21-24-55 


180 


741 


+89+55 
-89-49 


74fi +59+23+54 
'^"-52-23-53 


751 


+24+23+55 
-20-23-54 



and dileptonic decay channels, using 0.036fb . Both measurements are in agreement with 
our predictions within the uncertainties. 

5.3 Top-quark mass determination 

The ATLAS experiment also extracted the top mass from the measurement of the total 
cross section [9J by comparing the NNLOapp prediction of [7j to the cross section mea- 
sured using different top-mass values in the kinematical reconstruction of the events. With 
35pb~^ of data and an experimental cross section au = 186.3 pb±5%(stat) ± 12%(syst) ± 
3%(lumi) for the reference top-mass rrit = 172.5 GeV, the ATLAS collaboration obtained 
the value rrit = I66.4I7 3 GeV for the pole mass. With a similar procedure, the DO collabo- 
ration obtained the pole mass [8] mt = 167.5^47 GeV from the cross section measurement 
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at the Tevatron. 

It is interesting to speculate how the central value and the accuracy of the mass determi- 
nation will change with the improved statistical accuracy of the latest LHC measurements 
and our NNLL2 predictions. Following the method described in |9j we define a likelihood 
function 

fim) = J fthio-\mt) ■ fe^p{a\mt)dcr , (5.2) 

where /th is a normalized Gaussian distribution centred on the theoretical prediction of 
the cross section, having a width equal to the total error of the theory prediction 



/th(o-|mt) 



1 



exp 



2(A<(m,))2 



(5.3) 



/exp is defined in a similar way, with the central value and width of the Gaussian given 
by the measured value of the cross section, a^^^{mt), and the total experimental error, 
A(T^f^(mt). The top-quark mass is then extracted from the maximum of the likelihood 
function (15. 2p . with the error obtained from the 68% area around the maximum. 

The theoretical cross section ct*j (mj) is obtained using the numbers for the NNLL2 
approximation given in Table [TTl interpolated by a function of the form 



(172 5 \ ^ 
(co + ci(mi - 172.5) + C2{mt - 172.5)^ + c^{mt - 172.5)^) pb , 

(5.4) 

where all the masses are given in GeV. The total theory error is obtained by summing in 
quadrature the scale, resummation and the NNLO-constant uncertainties. For the PDF+a^ 
uncertainty, the 68% CL set was used in the Atlas note |9] and added linearly to the theory 
error. We account for this by rescaling our PDF-l-a^ error. Furthermore, when the errors 
in the positive and negative directions differ from each other, the maximum of the two 
values is used. For the coefficients in the fit (15. 4p we find cq = 166.496 ± 7.706 ± 7.464, 
ci = -(1.15093 ± 0.05687 ± 0.03741), C2 = (5.06265 ± 0.15147 ± 0.15886) x lO^^, C3 = 
(8.53722 ±12.1116 ±(-3.19098)) x 10"^ where the first error denotes the total theory error 
and the second the 68% CL PDF±as error. Using the theoretical prediction provided by 
our NNLL2 approximation, the measured ti cross section for the 35pb~^ data set given in 
Table 1 of |9], and the extraction procedure described above, we obtain mt = 166.5l7;g GeV, 
which is in good agreement with the value extracted by the ATLAS collaboration using 
the NNLOapp result from [7]. 

For the more recent measurement obtained by ATLAS with improved statistics, the 
dependence of the measured cross section on the reference top mass was found to be well 
described by a linear fit ^ 



a^^'^imt) = (411.9 - 1.35 mt) pb 



(5.5) 



in the interval rrit = 160 — 190 GeV. The mass dependence of both experimental and 
theoretical cross section is shown in Figure [61 Using the experimental result (15. 5 p and our 
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o-tt(m,)[pb] 




Figure 6: Mass dependence of the theoretical NNLL2 cross section (red) and of the measured 
cross section (black), as obtained from Ref. [5]. The solid lines represent the central values, 
while the total uncertainties of the theoretical and experimental results, determined as explained 
in the text, are given by the external dashed lines. 

theoretical prediction (15 .4^ for the evaluation of the likelihood function, we extract the 
pole mass 

= (169.8tt?)GeV, (5.6) 

where we have assumed that the total relative experimental error does not depend on m^, 
and is equal to the relative error at the reference mass, Aa^^^/a^^^ = ±6.6%. It can 
be seen that the central value obtained from the 0.7 fb~^ data set is higher than the one 
obtained in [9], and the error is reduced to ±3%. Also note that (15. 6p agrees with the 
Tevatron measurement from direct reconstruction, rrit = 173.3 ±1.1 GeV, at better than 
la accuracy. 

6 Conclusion 

We calculated the total production cross section of top-antitop pairs at the Tevatron and 
the LHC based on a combined resummation of threshold logarithms and Coulomb correc- 
tions at NNLL accuracy, using the formalism developed in [3ll|39]. We carefully assessed 
the ambiguities inherent in the resummation prescription, and adopted a new procedure 
for choosing the soft scale in the momentum-space approach to soft-gluon resummation 
in addition to the fixed scale advocated in (SUES], obtaining consistent results with both 
methods within the estimated uncertainties. We also used the expansion of the NNLL 
result to obtain new approximate results at N^LO accuracy. 

Our calculation accounts for bound-state corrections and higher-order Coulomb cor- 
rections not included in NNLL or NNLOapp calculations based on summing threshold 



37 



logarithms at fixed invariant mass [STJ [3811111112]. We find that the main effect of these 
corrections is included in the NNLOapp corrections obtained in [36], whereas the Coulomb 
corrections beyond NNLO are in the sub-percent range. The effect of NNLL resummation 
compared to the NNLOapp result of [36j is small, about 2% at the Tevatron and below 1% 
at the LHC While this suggests that effects beyond NNLO are small, it is interesting to 
note that a larger effect is obtained by adding the partial N^LO results obtained from the 
expansion of the NNLL prediction. This might indicate that, from a certain perturbative 
order onwards, resummation should always be performed, even though the total effect is 
small. 

Our results for a pole mass of nit = 173.3 GeV, and using the MSTW08NNLO PDFs, 
are given by 



where the first error denotes the combined theoretical uncertainty, including the scale 
uncertainty, estimates of resummation ambiguities and the unknown constant term in the 
threshold expansion of the NNLO cross section, and the second the PDF+as uncertainty. 
The first error should be reduced once the full NNLO calculation of the top cross section 
is available. The reduction of the second one depends on the availability of a suitable 
reference cross section that depends on the gluon luminosity. 

Using our NNLL2 prediction and recent ATLAS measurements, we estimated that the 
top-quark mass could be extracted with an accuracy of ±5 GeV from the currently available 
LHC data on the total cross section, and our result 



is compatible with the mass determination from direct reconstruction. 
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(Ttf(Tevatron) 
a«-(LHC, : 
a«-(LHC, : 




7TeV) = (162.61^:^™) pb 
14TeV) = (896^^°^^^) pb 



(6.1) 



mt = (169.8^1?) GeV 



(6.2) 



38 



A NLO Coulomb function 



The Coulomb Green functions at leading and next-to-leading order appearing in (12.1 ip 
read ESl''^ 



= 'it 



m 



t_ 
An 



rrit 



'-DR)as 



L^- 1+^(1 -A) 



4V) = 1^4vr(-D,K^ 



ai [Le + Jo] + /3o [Ll + '^3()Le + Ji] | , 



(A.l) 



where we have introduced the modified Euler psi function, ip{x) = 7^ + '^/'(a;), the short- 
hands 

Lk = --In 



A 



-E/rut 



and the functions 

Jo = AV^'(1 - A) - - A) 



ji = 44F3(1,1,1,1;2,2,1-A;1) + A^"(1-A)-2A^(1-A)V''(1-A) 



-3^'(1-A) + ^(1-A)'- — 



6 



(A.2) 



(A.3) 



(A.4) 



The coefficients /So = 11 — -^rii and ai = ^ ^ni {rii = 5 light quarks) denote the first 

coefficient of the QCD beta-function and the one-loop correction to the Coulomb potential, 
respectively. The computation of the convolution fl2.18p requires the numerical evaluation 
of the function ji and its derivatives for complex argument A, which is explained in Section 
lA.ll After the analytic continuation of the convolution of partonic cross sections with 
the parton luminosity, the potential function from f l2.1ip with set to one needs to be 
evaluated at = 0, see (I2.19p . For the colour- singlet configuration {Dr = —Cp) this is 
given by 



Ji(0) 



2ti 



a. 



(for A 



nC 



1). 

(A.5) 



The octet Coulomb function vanishes in this limit. 



A.l Evaluation of the NLO Coulomb function 

The only non-trivial term in the NLO Coulomb function is the hypergeometric function 
4F3 in (lA.4p . which needs to be evaluated for complex A. We define 



F43= 4^3(1, 1,1,1; 2, 2, JV;i; 



E 

i=0 



T{t + lfT{N) 



with = 1 - A , 



(A.6) 



^The expression for the NLO Coulomb- Green funetion obtained in [53] is quoted e.g. in [56] . 
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By assuming that is a positive integer, this sum is of the type which has been considered 
in [63], and can be expressed in terms of harmonic sums [781179] . 



(A.7) 



F43 = (N-l) jCa - C2Si{N - 2) + Ss{N - 2) + S,{N - 2)S2{N - 2) 
-S2,i(iV-2)| . 



The single harmonic sums can be continued analytically to complex by expressing them 
in terms of the ■j/'-function and its derivatives 



Si{N) = ^(Ar + l)+7^ 



Sa{N) 



r a) 



(a-l) 



(A^ + 1) + Ca , for a G N, a > 2 . 



(A.8) 
(A.9) 



The continuation of the nested harmonic sum 5*2,1 (A^) is well-defined as well and can be 
easily implemented, cf. e.g. [80l|8T]. The result for F43 in terms of A G C then reads 



^43 



C2 - 52(-A) - A Cs + S^{-X) - 5i(-A) (C2 - S^i-X)) - 52,i(-A) 



(A.IO) 



The computation of the convolution (12.181) requires the numerical evaluation of derivatives 
of the NLO Coulomb function with respect to A. Derivatives of the ?/^-function and hence 
single harmonic sums are easily done, but the derivative of the S'2,1 requires more work. We 
use the package HarmonicSums, [82] for this task and obtain for the first derivative of F43 

^F43 = -a| + C2^i(-A) + AC252(-A)-^i(-A)52(-A) 

-\^^^ - 3S3(-A) - 2XS^{-\)S,{-\) - 5\^^ 

+^2,i(-A) + 2A^3,i(-A) + C3 + 2A^i(-A)C3 . (A.ll) 

In ( lA.lip an additional nested harmonic sum appears, 5*3,1, which can again be continued 
to complex A using methods described in [SOllSTj . 



B Bound-state corrections 

Above the production threshold, E > 0, the potential function J/j is determined by the 
branch cut of the Coulomb Green function Gc\r{E) along the real axis. Since we neglect 
the top decay width in the potential function, which is allowed when computing the total 
cross section, an imaginary part arises below threshold from isolated poles corresponding 
to ti bound-state production. The bound-state contribution to the potential function takes 
the form 

00 

jR{E)=2j2HE-E^)R^e{-Dn), E<(} (B.l) 

n=l 
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There are no bound states corrections when the Coulomb potential is repulsive {Dn > 0), 
as is the case in the colour-octet production channel, as indicated by the step function 
9{—Dfi). The leading-order bound-state energies and the residues read 



^(0) _ aiDw^t 



mt{-DR^)as 
2n 



(B.2) 



While these expressions are sufficient for NLL resummation, the NLO Coulomb Green 
function shifts the energy and residue of the bound states to 



with EaiE 



ei = 2ai + 4/3o 



6ri = 3ai + 2(3q 



Siiji) — In 



rritas 



-D, 



Siin) + 2nS2in) 



nrr 



3 In 



rritas 



(B.3) 



(B.4) 



(B.5) 



The bound-state contributions to the partonic cross section are convoluted with the 
parton luminosity, similarly to the procedure discussed in Section 12.31 We mention here 
the relevant modifications. Recall that we do no apply the small correction Anc in (12. lip to 
the already small bound-state contribution. Inserting fIB.ip into the resummation formula 
and the convolution with the parton luminosity, we find 



a 



BS _ 
NiN2^ttX — 

E 

p,p'=<i,g<9 

X 



J2 H^Amt,fi) / drLppir) / du^E - u/2)W\u, fi)e{u/2 - E) 
J ^„ Jo Jo 



p,p=q,q,g 
2Hl ,{muij) 



n 



rrit 
1 



exp[-4S'(/x^,/x/) + 2aJy(/x^,/i/)]s^(9^,/Xs) 



r(2r/) V ^^ 



2rj 



dTLpp,{T) 



E — E„ 
nit 



2»?-l 



(B.6) 



where E = y/rs — 2mt, and Hpp,{mt, fi), §■'■(9^, /is) are the hard function and Laplace 
transform of the soft function for the singlet state. The integrals in the series have different 
lower integration limits, given by 



1 2 

Tn = - {2mt + En) < To 



Am. 



(B.7) 



s s 
For r/ < the r integrals have a non-integrable singularity at r = r„ where E — En 



ys/(4r„) (r — Tn). This is analogous to the convolution of the continuum cross section 
with the parton luminosity in Section 12.31 Following the method discussed there, the 
analytic continuation to negative r] is constructed. However, here the integrand behaves 
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~ (t — TnY^^^ rather than ~ (r — tq)^'', so one further subtraction is required. To improve 
numerical stabihty we again introduce a parameter A > tq and spht the integral into two 
integration domains, 



/"i / — E ^ 



nit 

2-q-l 



^-^ \ nit 



+ £dTLppiT) (^^^ ' 5^i?„ {^-^f-' . (B.8) 



In the second term we interchanged summation and integration. The second integral is 
free of singularities. Furthermore, the sum inside the integral is convergent. To isolate 
the singularities in the region [r„. A] we introduce the expansion of the integrand around 
r = Tr, as 



,2r,-l 



] = (47 



A/2-ri 



LppiTr,)iT-T^Y^-' 



ATnL'^,{Tn) + (1 " 2r])Lpp>{Tn) 



4r„. 



[T - Tr, 



(B.9) 



Then, for i] < 0, the integrals over the interval [r„. A] are analytically continued to 77 > — 1 
by the identity 



[ drLppir) {VT - ^Y'-' = [ dT{l-r}!^)[LppiT){VT 



Lpp'{Tn) , 4r„L; + (1 - 2T])Lpp,{Tn) ' 

[A-Tn) 



2r] 



4(1 + 2r])Tn 



. (B.IO) 



The poles at 1] = 0, —1/2 are now manifest and cancelled by the factor l/r(2?7) in flB.6p . 
The remaining integral in fIB.lOp is singularity- free if 77 > —1, and can be computed 
numerically. This requires the numerical evaluation of the first derivative of the parton- 
luminosity functions, L'pp,. The integral fIB.lOp can be continued to arbitrary negative 
values of i] by performing more subtractions. However, for the applications presented in 
this paper, the continuation to rj > —1 is sufficient. The convergence of the series in n, and 
the stability of the numerical integration, depend on the choice of the separation parameter 
A. We adopt A = 1.002ro, and truncate the series over bound states at level nmax = 100. 



C Fixed-order expansions 
C.l Expansion to 0{a1) 

In this appendix we provide the NNLO corrections to the cross section (13. ip resulting from 
expanding the NNLL resummed result. The scale- independent terms for the three partonic 
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channels are given by 

4-(8) = " 02 + « (-140.368 In^ /3 + 32.106 In /3 + 3.95105) 

+ 910.222 In^ (5 - 1315.53 W ^9 + 592.292 In^ (3 + 528.557 In /3 + /^^Jgj , (C.l) 



/£(!) = ^^^^ + -R (2526.62 In^ /3 + 1347.76 In /3 + 66.1114) 



230.896 1 

+ 4608. In^ /3 - 249.2 In^ /3 - 2385.73 In^ /3 + 1600.47 In /3 + /^^^l} , (C.2) 
/^2,o) _ • + (-315.827 In^ /3 - 89.5134 In /3 - 38.5162) 

+ 4608 In^ (3 - 2553.2 In^ /3 - 322.996 In^ /3 + 2799.24 ln/3 + /^^^^J^ ■ (^-3) 

The terms shown exphcitly are known exactly and have first been obtained fully in |36] . 
The functions /^'*^^ are constant terms in the threshold expansion that are currently not 
known completely, since, for example, the two-loop hard functions contribute to them. At 
the NNLL level, they therefore show a residual dependence on the numerical factors ki in 
the scale choices /i^ = ksmtf3^, and fih = khmt- These terms are given by 

= - 56.889 In^ kh + 240.692 In^ kh - 116.423 In^ kh - 456.803 In kh 

- 56.889 In^ k, - 164.441 In^ k^ - 52.0528 W k^ + 319.029 In k, 

+ 49.7741 , (C.4) 

•^^0)^^^'°^ = - 288. In^ kh + 921.172 In^ kh - 8.19545 In^ kh - 1322.62 In kh - 288. In^ ks 
31.15 In^ k, + 1231. In^ k^ + 1295.98 In ks - 328.235 , (C.5) 



^NLL(2,o) ^ _ 2gg_ _^ 633.172 In^ kh + 543.378 In^ kh - 1109.7 In fc/, - 288. In^ fc, 

- 319.15 In^ ks + 847.638 W fc, + 1626.13 In ks - 43.6962 . (C.6) 

The factorization-scale dependent terms for the partonic channels are 

/^2>i) ^ 1 (403.557 - 2526.62 In (3) - 9216 In^ /3 + 3567.82 In^ /3 

+ 2639.47 In + J), (C.7) 

=46081n2/3 - 3563.961n/3 + /^^J), (C.8) 
= 1 (315.827 In /3 - 50.4446) - 9216 In^ /3 + 5871.82 In^/? 

-25.28771n/3 + 4Y8)' (^.9) 
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In this case the constant terms are known exactly: 

4-(8) = 327.048, /(Jg = 385.383, (C.ll) 

= -487.213, /(^g = -417.165, (C.12) 

& = -283.914, fllf^, = f!^lfl (C.13) 

C.2 Expansion to 0{al) 

Here we provide those parts of the N^LO corrections to the cross section (13. ip that have 
not been given aheady in Section [31 i.e. the scale dependent terms f^'^^ for i = 1, 2, 3 and 
the NNLL approximations to the functions fx'^\ that are beyond the NNLL accuracy of 
the threshold expansion. The expressions that follow are generated from expanding the 
resummation formula as given in Eqs. (3.48) and (3.49) of [39] (appropriately adapted to 
the case of top-quark production). 



Quark- antiquark channel. The scaling functions of the scale-dependent contributions in 
the N^LO expansion of the quark-antiquark production channel are given by 



f(3,l) 
Jqq{8) 


1 


r(-153.931n/3 + 56.8546) + ^ (5989.02 In^ /3 - 7733.19 In^/? 


- 2669.2 In /3) 






38836.1 In^ /3 + 109310. In^ /3 - 78403.7 In^ /3 + /^^Jgj , 


(C.14) 


f(3,2) 
Jqq{8) 


1 


(-2994.51 In^ /3 + 5287. 18 In /3) + 38836.1 In^ /3 - 111164 In^ /3 + /^^Jgj , (C.15) 


f(3,3) 
Jqq{8) 




12945.4 In^/? + 4^3;, 


(C.16) 



where all the terms shown explicitly are exact. These are the only terms included in the 
N^LOb approximation, while the N^LOa approximation includes in addition the NNLL 
approximation to the /(^'*)-functions. Most of the terms in these functions can only be 
predicted exactly from resummation beyond NNLL accuracy. At the NNLL level, they 
show a residual dependence on the numerical factors kt in the scale choices fig = ksmt/3'^, 
fih = khirit and /ic = kcrritP. For the quark-antiquark channel these functions read 

/^^^^'"^ = ^ (773.485 In^ fee - 335.834 In fc^ + 187.157 In^ 4 - 791.846 In^ fc;, 

+383.015 In^ A;/, + 1502.82 In A;,, + 187.157 In^ A;, + 81.5526 In^ fc, 
-1199.28 In^ k, - 2174.59 In k, + 749.55) + In /3 (-2331.66 In kc 
+1172.51 In^ kh - 4960.79 In^ kh + 2399.53 In^ kh + 9414.94 In kh 
+4661.69 In^ ks + 4398.91 In^ ks - 27481 In^ k^ - 60914.9 In ks + 11370.2) 
+ In^ (3 (-2427.26 In^ kh + 10269.5 In^ kh - 4967.37 W kh 
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-19490.3 In kh - 2427.26 In^ ks - 7016.16 In^ ks - 2220.92 In^ ks 
+13611.9 In ks - 65560) - 404.543 In^ kh + 3730.45 In^ kh 

- 10015.2 In^ kh + 5148.7 In^ kh + 9150.36 In^ kh - 5460.36 In kh 
+ 404.543 In^ ks + 2917.1 In^ ks + 1022.1 In^ ks - 12319 In^ ks 

- 15042. lln2A;s + 5890.83 In A;, + 702.307, (C.17) 

f^m^^^'^^ ^ 22^:92 _ 3gQ^2.7 In^ /3 + In /3 (2427.26 W kh - 10269.5 In^ kh 

+4967.37 In^ kh + 19490.3 In kh + 2427.26 In^ A;, + 7016.16 In^ ks 

+2220.92 In^ A;, - 13611.9 In A;, + 67825.2) - 3323.77 In^ A;^ + 14062.6 In^ A;,, 

- 6802.08 In^ kh - 26436.3 In kh - 3323.77 In^ A;, - 5069.58 In^ A;, 

+ 12832 In^ A;, + 36025.7 In A;, + 327.661 , (C.18) 

^T^NLL(3,2) ^ _ 1005-67 ^ ^g^^Q ^ _ 3333 32 In ^ - 3697.7 , (C.19) 

fgm^^^'^^ = 39223.6 W /3 - 28867.1 \n/3 + 1140.45 . (C.20) 

Gluon fusion colour- singlet channel. The scaling functions for the factorization-scale de- 
pendent terms are 

/^(^>i^) = (-22166 In ^ - 8283.49) + ^ (-242555 In^ ^ - 51902.1 In^ ^ + 217024 In 0) 

- 442368 In^ /3 + 300977 In^ /3 + 563908 In^* /3 + /^^^ , (C.21) 

/^(3:2) _ ^ (121278 In^ ^ - 58112.2 In + 442368 W ^ - 448309 ln=^ /3 + /^^^^^^^ ' (^-22) 

/Si) = - 147456 W ^ + . (C.23) 
The NNLL approximations for the terms in the scaling functions not exactly known read 

7^(3:0) _ 1 (-6187.88 In^ kc + 2686.67 In A;^ - 7579.86 W kh + 24244.3 In^ kh 

-215.695 In^ kh - 34809.9 In kh - 7579.86 In^ ks + 17787.4 In^ ks 
+49075.2 In^ ks + 10292.7 In ks - 14467) + In /3 (36589.2 In kc 
-2196 A In^ A;/, + 7025.22 In^ kh - 62.5017 W A;,, - 10086.8 In kh 
+15467.6 In^A;^ - 101709 In^ ks - 211987 In^ ks - 10792.4 In A;, + 112623.) 
+ In^ (3 (-27648 In^ kh + 88432.5 In^ kh - 786.763 In^ kh - 126971 In kh 
-27648 In^ A;, - 2990.4 In^ ks + 118176 In^ ks + 124414 In A;, - 146090) 

- 4608 In^ kh + 27996.1 In^ kh - 42685.9 In^ kh - 17071.1 In^ kh 
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+ 45365.4 In^ kh + 13255 In kh + 4608 In^ ks + 6635.6 In^ ks - 62517.1 W ks 



79096.6 In^ ks + 37998.9 W ks + 37910.9 In ks - 9739.8 , 



(C.24) 



?-NNLL(3,l) 
J 99(1) 



29709.3 



- 618319 In^ /3 + \n/3 (27648 In^ kh - 88432.5 In^ kh 



NNLL(3,2) 



+786.763 In^ kh + 126971 In kh + 27648 In^ ks + 2990.4 W A;, 

-118176 In^ A;, - 124414 In ks + 198168) - 17315.9 In^ kh + 55385.1 In^ kh 

- mM7 In^ kh - 78953.1 In kh - 17315.9 In^ ks + 49817.9 In^ A;, 

+ 146950 In^ ks + 48514.7 In ks - 77934 , (C.25) 

43685.6 



297750 In^ p + 322713 In ^ - 52925.9 , 



f 99^(1)^^^^^^ = 206398 In^ /3 + 10843.3 In /3 - 50289.8 . 



(C.26) 
(C.27) 



Gluon fusion colour-octet channel. The scaling functions for tiie factorization-scale de- 
pendent terms are 



1 1 
4.9(8) = ^ (-346.343 In /3 - 129.429) + - (30319.4 In^^ - 1092.09 In^ /3 - 21802.5 In /?) 

(C.28) 



'flfl(8) ^2 V ^ ' ^ 

- 442368 In^ ^ + 522161 In^ ^ + 227359 In^ /3 + /^^^gj , 
1 



ffnfl = 4 (-15159.7 In^ P + 7264.03 In /3) + 442368 In^ ^ - 558901 In^ /3 + /^^^|} , (C.29) 
«=- 147456 ln-^;3 + 4,). 



'ffff(8) " 

[.(3,3) _ _ 1 ^7/1 i„3 I f(3,3) (C 30) 

The NNLL approximations for the terms in the scaling functions not exactly known read 



•NNLL(3,0) 
99(8) 



?1SrNLL(3,l) 
■^99(8) 



^ (773.485 In^ kc - 335.834 In kc + 947.482 In^ kh - 2083.05 In^ kh 

-1787.64 In^ kh + 3650.77 In kh + 947.482 In^ ks - 1275.94 In^ ks 
-6617.62 In^ ks - 4868.71 In ks + 2424.05) + In^ /3 (-27648 In^ kh 
+60784.5 In^ + 52164.3 In^ kh - 106531 In kh - 27648 In^ ks 
-30638.4 In^ A;, + 81373.3 W ks + 156108 In ks - 347424) 
+ In /3 (-2152.3 In kc + 4715.6 In^ kh - 10367.3 In^ kh - 8897.06 In^ kh 
+18169.8 \\ikh + 22379.6 In^ ks - 78581.5 In^ A;^ - 289277 W A;^ - 205101 In A;^ 
+123426.) - 4608 In^ kh + 21084.1 In^ kh - 19739.3 In^ kh - 26067.8 In^ kh 
+ 33695.1 In^ kh + 1303.32 In kh + 4608 In^ A;, + 13547.6 In^ A;, - 50831.9 W A;^ 
- 118982 In^ A;, - 16519.8 In^ A;^ + 62214.4 In A;^ + 966.578 , (C.31) 

5903.02 



- 748536 W (5 + hi 13 (27648 In^ kh - 60784.5 In^ kh 
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Figure 7: Determination of /3cut and the resulting ambiguity for the Tevatron with kg = 2. The 
plot shows att{A^, i?>, /3cut) (in pb), as defined by convoluting (14.141) with the parton luminosity. 
Solid curves: = NNLL2, dashed curves A< = NNLLi, black curves: fi> = NNLL2, red 

N^LOa. The 



curves: 



B 



> 



NNLO 



app. 



blue curves: By = N'^LOb, green curves: B 



> 



horizontal dotted line is the NNLOapp cross section, the x-axis intersects the y axis at the value 
of the NLO cross section. The point with the error bar is the central value of the cross section 
together with the error estimate, obtained by taking the envelope of the curves in the ±20% 
window around the default /3cut = 0.35, given by the vertical dashed lines. 



'99(8) 

7NNLI 
'99(8) 



7NNLL(3,3) 
J n 



-52164.3 In^ h + 106531 In kh + 27648 In^ k^ + 30638.4 In^ ks 
-81373.3 In^ ks - 156108 In fc, + 399352) - 17315.9 In^ kt + 38069.2 In^ kh 
+ 32670.3 In^ kh - 66151.5 In 4 - 17315.9 In'* k^ + 32502.1 In^ ks 
+ 162669 W ks + 137140 In ks - 61943.1 , 



/NNLL(3,2) ^ 5^7 _ iQQ578i^2^ ^ 284537 In /3 - 88163.8 



= 1 



^NLL(3,3) 
99(1) 



(C.32) 

(C.33) 
(C.34) 



D Details of Method 2 

We provide here the numerical details defining /3cut in Method 2, as discussed in Section IT2| 
and the corresponding implementation at the NLL order. Figures [7] and [8] show the various 
approximations fl4.14p for /c^ = 2 at the Tevatron and at the LHC, respectively, together 
with the default /3cut and the 20%-window used in the estimate iii) of the resummation 
ambiguity. The plots show the expected behaviour: for increasing /3cut (and hence fif) 
resummation becomes less effective and the difference between the two implementations 
used in the lower interval (the difference between solid and dashed lines of the same colour) 
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o-(A,B,/3cut) 




Figure 8: Determination of /3cut and the resulting ambiguity for the LHC with 7 TeV (above) 
and 14 TeV (below). The curves are the cross section crtj(y4<, i?>, /3cut), with the various 
approximations defined as in Figure [71 

becomes dominant. For smaller P^nt the difference of the resummed cross section and 
perturbative expansions to different orders used in the upper interval (e.g. the difference 
among the solid lines of different colour) becomes sizeable, while the difference between 
the two NNLL implementations in the lower interval becomes negligible. Note that the 
default curve (the black solid hne) depends only weakly on the precise value used for 
/Scut- The spread among the different approximations for varying P^ut by 20% leads to an 
uncertainty of ~ 2.5% for both Tevatron and LHC. This is the dominant contribution to 
our estimate of the resummation uncertainty. The criteria i) and ii) discussed in Section 
both contribute uncertainties ~ 1.5% or smaller, cf. Table [3]. 

Finally we also present the values of Pcut one obtains when varying ks by a factor of 2 
around its default value kg = 2 for = 173.3 GeV. We obtain 
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/3cut(NNLL) : 0.52 (Tevatron), 0.75 (LHC7), 0.79(LHC14), 

/i< = : 47 GeV (Tevatron) , 97 GeV (LHC7) , 108 GeV (LHC14) . 



(D.l) 



/3cut(NNLL) : 0.23 (Tevatron), 0.33 (LHC7), 0.33(LHC14), 

jjf = AmtPl^^ : 37 GeV (Tevatron) , 76 GeV (LHC7) , 76 GeV (LHC14) . 



(D.2) 



As one would expect, the corresponding soft scales (if lie not very far apart from each 
other for fc, = 1, 2, 4, Eq. fHTT^ for k2 = 2. 

For resummation at NLL order, we determine Pent anew, using the same prescrip- 
tion to minimize the spread between different cross sections fl4.14p for the approxima- 
tions A< G {NLLi,NLL2} (defined in f lOj) and ( 144|) ) in the lower interval and i?> G 
{NLL2, NLO, NLL-NNLOa, NLL-NNLOb} in the upper interval. Here, in analogy to the 
two N'^LO approximations, NLL-NNLOa denotes the expansion of the NLL corrections to 
order a^, keeping all the generated terms, while only those terms exactly reproduced at 
NLL accuracy (i.e. the terms { j^, In^ /3, In^ /?}) are kept in the NLL-NNLOb 

approximation. With this prescription, the resulting /3cut-values at NLL for kg = 2 are 
somewhat smaller than for NNLL resummation: 

/3cut(NLL) : 0.26 (Tevatron), 0.32 (LHC7), 0.32 (LHC14) (D.3) 

The remaining uncertainties are estimated as for NNLL. 
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